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ABSTRACT 

We conduct an ultraviolet {HST and FUSE) spectroscopic survey of H I (Lyman lines) and seven 
metal ions (O VI, N V, C IV, C III, Si IV, Si III, Fe III) in the low-redshift intergalactic medium 
(IGM) at z < 0.4. We analyzed 650 Lya absorbers over redshift pathlength Az = 5.27, detecting 
numerous absorbers: 83 O VI systems, 39 C III, 53 Si III, 24 C IV, 24 N V, and so on. Our survey 
yields distributions in column density and estimates of the IGM baryon content and metallicities of 
C, N, O in the IGM. In the low-z IGM, we have accounted for ~ 40% of the baryons: 30% in the 
photoionized Lya forest and 10% in the (T = 10^~^ K) warm-hot intergalactic medium (WHIM) traced 
by O VI. Statistical metallicities of C, N, O ions are consistent with the canonical (z = 0) value of 10% 
solar, with considerable scatter. Improved statistics for weak O VI absorbers allows us to estimate 
^wmm/^b = 0.073 ± 0.008 down to log A^ovi = 13.4 and 0.086 ± 0.008 down to log Nqyi = 13.0. The 
O VI absorber line frequency, dM jdz = 40lg"', down to 10 mA equivalent width suggests a 250-300 
kpc extent of metals around dwarf galaxies (My = —17 to —18). Many absorbers appear to contain 
multiphase gas, with both coUisional ionization and photoionization determining the ionization state. 
N V absorption is well-correlated with O VI and both ions show similarly steep power-law indices 
dM /dz (X N^f^ with /3ovi ~ /3nv ~ 2 while /3hi = 1-7. We conclude that O VI and N V are reliable 
tracers of the portion of the WHIM at T « 10^~^ K. C IV may be present in both coUisional and 
photoionized phases; A^civ correlates poorly with both A'hi and A"ovi and /3hi < /3civ < /3ovi- The 
ions C III, Si III, and Si IV are well correlated with H I and show patterns typical of photoionization. 
We detect Fe III in 14 absorbers with column densities higher than expected, [Fe Ill/Si III] 0.5±0.3. 
Adjacent ion stages of the same element (CIII/IV and Silll/IV) provide useful constraints on the 
photoionization parameter, \ogU « —1.5 ± 0.5. Comparison of Si IV and C IV with high-z surveys 
shows a modest increase in line density, consistent with increasing IGM metallicity at recent epochs. 
Subject headings: cosmological parameters — cosmology: observations — intergalactic medium — 
quasars: absorption lines — ultraviolet: general 



1. INTRODUCTION 

We continue our far-ultraviolet (FUV) surveys of the 
baryon content and metallicity of the low-redshift in- 
tergalactic medium (IGM) by analyzing 650 Lya ab- 
sorbers, the largest database to date from the Space Tele- 
scope Imaging Spectrograph (STIS) aboard the Hubble 
Space Telescope {HST) and the Far Ultraviolet Spectro- 
scopic Explorer (FUSE). These spectrographs offer sensi- 
tive probes of H I (Lyman lines) together with associated 
metal-line absorbers. Of particular importance are UV 
absorption lines from high ion stages (O VI, N V) pro- 
duced primarily in hot, coUisionally ionized gas, as well as 
lower ion stages (C IV, C III, Si IV, Si III, Fe III) that can 
arise in both photoionized gas and hot gas. The advent 
of these sensitive FUV spectrographs has stimulated an 
increase in our understanding of the low-redshift IGM. 

In this survey, we address several astronomical ques- 
tions facilitated by a more accurate census of baryons 
and heavy elements in the IGM. From obser vations with 
HST FUSE and Chandra, it seem s likely (" Skllil [20031 : 
IStocke et ani2006at iNicastro etldl [2005) that most of 
the baryons still reside in the IGM, even at low red- 
shift. But in what thermal phase and ionization state 
do they exist, and with what spatial distribution? What 
are the physical properties (temperature, density) of the 
low-z absorbers? From their metallicities and nuclcosyn- 
thetic patterns, can we uncover their stellar sources? By 



comparing the spatial proximity of metal-line absorbers 
to neighboring galaxies, can we estimate the extent of 
metal transport into the IGM? In this paper, we attempt 
to make progress on all of these issues. 

Over the p ast decade, the Colorado group and 
others (e.g., iTripp et"al1 l2008f) have focused on 

UV studies of H I and O VI absorb ers in 

the low-z IGM (iPenton. Stocke. fc Sl^ [Mil [20M 
iTripp. Savage, fc Jenkins! [20001: [Savage et al.l [2002f) . 
Danforth & Shull (2005) and Danforth et al. (2006), here- 
after denoted Papers 1 and 2, respectively, analyzed IGM 
absorption at z < 0.3 in FUSE sight lines toward 31 AGN 
with high-quality data from STIS and FUSE. Those stud- 
ies began with a set of 171 strong Lya absorbers (equiv- 
alent widths Wx > 80 mA) taken from the literature 
([Penton. Stocke. fc Shul l 2000, 2004, and other sources). 
Because Lya lines are subject to instrumental broad- 
ening and line saturation, we endeavored to measure 
higher-order Lyman transitions (Ly/J-Lyr;) in the FUSE 
data at Zabs < 0.3 in order to obtain accurate curve-of- 
growth (CoG) solutions to the doppler parameter, 6hi, 
and neutral hydrogen column density, -/Vhi ■ As fi rst noted 
in FC/5'£; studies of Ly/3/Lya ([Shull et al.[[2000f) and con- 
firmed in Paper 2, the Lya measurements alone will 
systematically underestimate iVui, a discrepancy that 
increases in stronger lines, an d overestimate 6hi- Pa- 
per 2 also confirmed the result of [Penton. Stocke. fc Shull 
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(|2004[ ) that the column-density distribution of H I ab- 
sorbers follows a power law, dN/dNm oc N^f , with 
Pm = 1.68 ± 0.11. The CoG-derived 5-values show that 
most of the H I absorbers arise in gas with T < 10^ K, 
incompatible with their existence in a hot (shock-heated) 
IGM phase. A small population of broad, shallow Lya 
absorbers m ay probe trace amounts of H I in hot, h ighly 
ionized gas (|Richter et al.ll2004l : [Lehner et al.ll2007f) . 

Along with the H I lines, Papers 1 and 2 measured 
corresponding absorption lines of the metal ions O VI 
AA1031.926, 1037.617 and C III A977.020. We found O VI 
to be the most accessible tracer of the warm-hot ion- 
ized medium (WHIM) , thought to arise in shock-heated 
gas a t T = 10^-6 K (iDave et al . 1999; Gen & Ostriker 
Il999| ). The FUV hues of O VI are currently the best 
tracer of the "cooler" portions of these shocks, reach- 
ing a peak ionization fraction, /ovi ~ 0.22 at temper- 
ature Tmax = l O^ **^ K in coUisional ionization equi- 
librium or GTE (jSutherland &: Dopitaj Jl993) and simi- 
lar va lues in non-equilibrium cooling ( Gnat fc Sternberd 
l2007f) . Owing to the relatively high solar oxygen abun- 
dance, (0/H)q « 5 X 10^**, and the substantial oscil- 
lator strength of the stronger line of the O VI dou- 
blet (1031.9261 A, / = 0.1325), the O VI ultravio- 
let transitions are detectable down to column densities 
A^ovi ~ 10"'^'^ cm~^, more than two orders of magnitude 
below the current sensitivity of the X-ray absorption lines 
of O VII (21.602 A, / = 0.696) and O VIII (18.969 A, 
/ = 0.416) studied with Chandra and XMM-Newton. Al- 
though the {ls-2p) X-ray transitions are crucial for de- 
tecting the hotter WHIM at T > 10^ K, their intrinsic 
line strengths, /A, are much weaker than the {2s~2p) 
O VI ultraviolet transitions and the detection limits and 
velocity resolution of current X-ray instruments make 
them much less sensitive. 

Papers 1 and 2 represented a significant statistical im- 
provement over previous studies of small nu r nbers of sight 
lines (e.g. iTripp. Savage, fc Jenkins! l2000t ISavage et al.l 
|2002[) with 40 O VI detections and Az 2.2 vs. ap- 
proximately 10 O VI absorbers over Az ~ 0.5. The de- 
tection statistics were sufficient for us to draw a num- 
ber of intriguing conclusions. We found that A^ovi was 
poorly correlated with iVui, and the power-law distri- 
bution of O VI column density, /3ovi = 2.2 ±0.1, was 
significantly steeper than for H I. These results implied 
a multiphase IGM with hot, shocked WHIM (traced by 
O VI) kinematically associated with cooler, photoion- 
ized H I. We inferred a total cosmological mass frac- 
tion in the WHIM of at least f^wHiM = (0.0022 ± 
0.0003) [/i7o(^o/0.1Zq) (/ovi/0.2)]-i 4.8 ± 0.9% of 
the total baryonic mass at z < 0.15. Unlike O VI, the 30 
G III detections showed reasonable column density cor- 
relation with H I, and the dM jdz distributions showed 
similar slopes, suggesting that G III is predominantly 
photoionized. 

In this survey (Paper 3), we expand the parameter 
space of the low-z IGM survey, using a much larger data 
set of 650 Lya absorbers, observed over a total redshift 
pathlength Az = 5.27 (compared to Az = 2.2 in Pa- 
pers 1 and 2). There were several limitations inherent in 
our initial survey. Our combined STIS/E140M+f 
spectral survey provides an expansion and improvement 
in three important ways: 



(1) In Paper 1, we limited our search to 
"strong" Lyg systems (Wt.vq - > 80 mA). However, 
iPenton. Stocke. fc Shulll (|2004l ) showed that the median 
Wi,ya = 68 mA, so this equivalent-width threshold 
eliminated more than half of the possible IGM systems 
from study. Furthermore, Paper 1 found little if any 
correlation between iVni and A^ovi- Indeed, at low col- 
umn densities, logA^jji < 13.4, we found many absorbers 
with A^ovi ~ A^Hi- Thus, measurable O VI absorption 
is possible even in weak H I absorbers. We relax this 
"Lya bias" here by measuring metal ion absorption in 
all detected Lya systems. 

(2) The combined STIS/E140M+FC/S'£;dataset covers 
a large wavelength range : 905-1710 A in aU cases and 
out to 1729 A in most of the datasets. We use the longer 
wavelengths available in the E140M data to measure H I, 
O VI, and G III at higher redshifts than were possible in 
the FUSE da,ta, alone (zqvi < 0.15). We now cover Lya 
out to z < 0.4 and measure higher-order Lyman lines 
in strong systems to even higher redshifts in some sight 
lines. This expanded redshift space more than doubles 
the O VI absorber pathlength. 

(3) We use the broad spectral coverage of STIS to mea- 
sure metal ions not found in FUSE data. Lines of N V 
and G IV can exist at temperatures nearly as high as 
O VI, and they may provide additional tracers of the 
cooler WHIM. We also measure lines of Si IV, Si III, and 
Fe III, which trace predominantly photoionized material. 
By observing adjacent ionization states of the same el- 
ements (GIII/IV, Silll/IV), we can make direct ioniza- 
tion measurements without the complication of relative 
abundance uncertainties. 

In § 2, we discuss our selection of AGN sight lines and the 
rationale behind our choice of spectral lines. In § 2.2 and 
§ 2.3, we discuss our technique for Lya line identification 
and the measurement and verification of metal ion lines. 
In § 2.4 and § 2.5, we detail our analytic procedure and 
methods for deriving quantities of cosmological interest 
such as dN/dz, /3ion, metallicity, and r^ion- Our basic 
detection results are presented in § 3.1. We compare new 
H I, O VI, and G III results to those of previous work and 
discuss the new ions (G IV, N V, Si III, Si IV, Fe III) in 
this context. In § 3.2, we compare observed line ratios 
with a set of GLOUDY models to derive typical IGM 
ionization parameters and temperatures. We discuss the 
implications of these results to cosmology in § 3.3. In § 4 
we summarize our conclusions. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Sight Line Selection and Data Reduction 

We have performed detailed FUSE+'STl'^ analyses on 
sight lines toward 35 bright AGN. The primary objec- 
tive of this census is to remove as much "rich sight line" 
bias as possible, while accumulating the longest avail- 
able pathlength through the low-redshift IGM. There- 
fore, we included all available AGN with both FUSE and 
STIS/E140M data of reasonable quahty. The combined 
dataset covers 905-1710 A at a resolution of 20 km s^^ 
or better, at a signal to noise per resolution element 
{S/N)^cs ?J 5 in most cases. This gives us a potential 
redshift pathlength out to z > 0.4 in Lya, Ly/3, O VI, 
G III, and most of the other species, but only out to 
z w 0.1 in G IV. Five of our sight hues (3G 273, Akn654, 
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TABLE 1 
FUSE+STIS IGM Sight Lines 



order Legendre polynomials. For each 10 A segment, we 
calculated the local S/N from the standard deviation of 
the normalized continuum data about the fitted contin- 



Sight Line 


R.A. 


Decl. 


2AGN 


STIS 


FUSE Tapers 




(J2000.0) 


(J2000.0) 




(ksec) 


(ksec) I & IP 



Mrk335 
HE 0226-4410 
PKS 0312-770 
PKS 0405-123 
HS 0624+6907 
PG 0953+414 
Ton 28 
3C 249.1 
PG 1116+215 
PG 1211+143 
PG 1216+069 
Mrk205 
3C 273 
Q 1230+0115 
PG 1259+593 
PKS 1302-102 
Mrk279 
NGC 5548 
Mrk 1383 
PG 1444+407 
Mrk 876 
30 351 
H 1821+643 
Mrk 509 
PHL1811 
PKS 2155-304 
Akn 564 
NGC 7469 



00 06 19.5 

02 28 15.2 

03 11 55.2 

04 07 48.4 
06 30 02.5 

09 56 52.4 

10 04 02.5 

11 04 13.7 

11 19 08.6 

12 14 17.7 
12 19 20.9 
12 21 44.0 
12 29 06.7 

12 30 50.0 

13 01 12.9 
13 05 33.0 

13 53 03.4 

14 17 59.5 
14 29 06.6 
14 46 45.9 

16 13 57.2 

17 04 41.4 

18 21 57.3 

20 44 09.7 

21 55 01.5 

21 58 52.0 

22 42 39.3 

23 03 15.6 



+20 12 10 
-40 57 16 
-76 51 51 
-12 11 37 
+69 05 04 
+41 15 22 
+28 55 35 
+76 58 58 
+21 19 18 
+14 03 13 
+06 38 38 
+75 18 38 
+02 03 09 
+01 15 23 
+59 02 07 
-10 33 19 
+69 18 30 
+25 08 12 
+01 17 06 
+40 35 06 
+65 43 10 
+60 44 31 
+64 20 36 
-10 43 25 
-09 22 25 
-30 13 32 
+29 43 31 
+08 52 26 



0.025785 

0.495 

0.223000 

0.572590 

0.370000 

0.234100 

0.329700 

0.311500 

0.176500 

0.080900 

0.331300 

0.070846 

0.158339 

0.117000 

0.477800 

0.278400 

0.030451 

0.017175 

0.086470 

0.267300 

0.129000 

0.371940 

0.297000 

0.034397 

0.190000 

0.116000 

0.024684 

0.016317 



17.1 
43.8 

8.4 
27.2 
62.0 

8.0 
33.0 
68.8 
26.5 
42.5 

5.8 
62.1 
18.7 

9.8 
95.8 

4.8 
54.6 
69.8 
10.5 
48.6 
29.2 
77.0 
50.9 

7.6 
33.9 
10.8 
10.3 
22.8 



97.0 
33.2 
5.5 
71.1 
112.3 
72.1 
11.2 
216.8 
77.0 
52.3 
12.0 
203.6 
42.3 
4.0 
668.3 
142.7 
228.5 
55.0 
63.5 
10.0 
46.0 
141.9 
132.3 
62.3 
75.0 
123.2 
60.9 
44.3 



Y 

^ 2.2. Identification of Lya Absorbers 

We identified the associated higher Lyman lines and 
mcKil absorption systems using the Lya lines as "sign- 
post" . In the software, we displayed the normal- 
izedJ STIS/E140M spectrum in 10 A segments, mark- 
ing Ythe ISM lines and IGM line positions correspond- 
ing^ previously-identified Lya systems (e.g.. Figures [T] 
an4|2]). Any unidentified line was centroided, measured, 
andy" initially assumed to be Lya. For weak lines, in- 
cluffing broad shallow absorbers, we expect the optical 
depths to be in the ratio, T(Lya)/T(Ly/3) « 6.2. If the 
wa'^length and equivalent width at the position of the 
corf^sponding Ly/3 and higher-order Lyman lines were 
consistent with the Lya identification, we accepted the 
absorber and entered it into a database. For strong ab- 
sorlpors (Wlyq ^ 150 mA), we would expect to detect 
Ly/y or higher-order Lyman lines. Quite a few poten- 
tial^trong lines, initially assumed to be Lya, were ruled 
outyby this criterion. Line identification was an itera- 
tive^process: we assumed weak, unidentified lines to be 
Lyt^ unt il proven otherwise. More careful analysis occa- 



Sight line included in Papers I & II? 



Mrk 509, Mrk 205, and PG 1116-^215) are lacking the 
longest wavelength order of the STIS echelle data and 
truncate at A = 1709 A. The other 30 sight lines include 
the zeroth-order data at 1711 < A < 1729 A, which we 
include where possible. 

Of the 35 sight lines, we excluded six tar- 
gets (NGC 3516, NGC 3783, NGC 4051, NGC 4151, 
NGC 4395, and NGC 4593) at .zagn < 0.01 because 
the available redshift pathlength is minimal and any 
absorbers could be intrinsic to either our Galaxy or 
the AGN. We also excluded the high-redshift object 
HS 1700-1-6416 (z — 2.72) since its interpretation was 
greatly complicated by high- z EUV lines from s everal 
partial Lyman-limit systems (|Fechner et al . 1 120061 ) . Our 
final sample of 28 AGN sight lines is presented in Table 1. 
They cover a total Lya redshift pathlength of Az = 5.27 
with (zLya) ~ 0.14 and (zagn) = 0.21. 

We retrieved the FUSE data from the Multimission 
Archive at the Space Telescope Science Institute (MAST) 
archive and reduced them as in Paper 2. We binned the 
data by three instrumental pixels. The FC/S'i? instrumen- 
tal resolution of ^ 20 km s^^ typically represents 8-10 
pixels or ~ 3 bins. We used primaril y FUSE data from 
the highest-throughput LiFl channel (jMoos "eraI1[200l . 
but went to LiF2 or SiC channels when LiFl did not cover 
the region of interest or to provide a check on marginal 
LiFl detections. The archival STIS/E140M data were re- 
duced at the University of Colorado by S. V. Penton. We 
smoothed the data over three pixels to better match the 
nominal instrumental resolution of ~ 7 km s^^. Data 
from both instruments were normalized in overlapping 
10 A segments. We selected line-free continuum regions 
interactively and fitted the spectral continuum using low- 



sionally showed that lines identified as weak Lya were 
better explained as metal lines of stronger IGM systems, 
high-order Lyman lines of a strong higher-redshift sys- 
tem, intrinsic AGN features, high velocity components 
of Galactic lines, or instrumental artifacts. 

Since the rest-frame line density of FUV absorbers in- 
creases dramatically at shorter wavelengths, we began 
our Lya search at long wavelengths and worked our way 
to the blue. The identification process was carried out 
iteratively until all absorption features were either iden- 
tified or ruled out as not being Lya systems. In several 
cases, Lya lines were blended or redshifted beyond the 
long wavelength end of the data (1710 or 1729 A, Zabs > 
0.406 or 0.422, depending on sight line), and strong H I 
absorbers were then identified through their Ly/? absorp- 
tion. Our Lya identifications were informed by those 
of previous work, in particular iPenton. Stocke. fc Shul] 
(2000, 2004), but we searched each sight line and verified 
each absorber independently. 

Although some fraction of strong IGM systems are 
composed of multiple velocity components, this is not 
always apparent in the often-saturated Lya line profiles. 
In cases where higher Lyman lines showed multiple, un- 
ambiguously separable components in Lya or Ly/3, we 
treated each component as a separate absorber. How- 
ever, most profiles were consistent with single compo- 
nents or showed some asymmetry which we could not 
split unambiguously into multiple features. These cases 
(the majority) are treated as single absorbers. Our 
counting statistics are intermediate b etween the "compo - 
nent" and "system" nomenclature of [Tripp et all (|2008^ ■ 
but generally closer to the latter classification. 

In total, we identified 660 H I absorbers via Lya 
and/or Ly/3 absorption and found several discrepancies 
with published line lists. We truncated the line list at 
z = 0.4 since data quality falls rapidly beyond 1700 A. 
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Fig. 1. — IGM absorption in the z = 0.06808 absorber toward PG 0953+414 showing typical, normalized FUSE and 
STIS/E140M data. Lya and Ly/3 (top left) show strong, consistent profiles {Wi^ya = 284 ± 13 mA and Wi^yfi = 128 ± 16 
mA). GUI is detected but Si III and Felll are non-detections. High ions are depicted in the right hand panels: O VI A1032 and 
both NV lines (center left) are detected with consistent profiles. The O VI A 1038 line is blended with a strong H2 transition 
and is not shown. C IV shows noisy but consistent detections in both bands of the doublet as well. The two Si IV transitions are 
not shown. Each panel is centered at the redshifted wavelength of the transition and covers ±500 km s~^ in either direction. 
Other detected features in the data are identified. The "g:" prefix denotes a Galactic absorption line, while a numerical suffix 
denotes the redshift of an IGM line. The source channel (e.g. STIS/E140M or FUSE LiF2a) is indicated in the lower right. 



This left us with a total of 650 IGM systems, presented^ 
in Tabled 

2.3. Line Measurements and Consistency Checks 

The next step was to perform detailed measurements 
of 11 UV resonance absorption transitions of seven ion 

^ Truncated sample tables are presented here. Full versions can 
be seen in the electronic version. 



species (Table |3l), spanning a range of ionization lev- 
els and strengths. In particular, O VI and C IV are 
normally interpreted as tracers of hot gas, with peak 
ionization fractions (/ion) at temperature Tmax in col- 
lisional ionization eq uilibrium (CIE). As computed by 
Suth erland k Dopital (|T993), O VI has maximum CIE 
abundance of /ovi — 0.22 at logTmax = 5.45, N V 
reaches /nv — 0.24 at logT^ax = 5.25, and few = 
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Fig. 2. — Same as Fig. [l] for two absorbers in the Mrk876 sight line. The vertical dashed lines in the panels show the 
redshifted locations of the z = 0.11487 (left) and z = 0.11586 (right) absorbers. Both absorbers show weak, broad Lya profiles 
{Whya = 87 ± 13 and 74 ± 5 mA, respectively), and the expected Ly/3 equivalent widths are below the level of the noise and are 
not shown. CIII shows a clean non-detection at both redshifts. Data from the stronger SilV transition (A1394) is consistent 
with noise and a better upper limit is given by the weaker SilV A1403 line. Si III shows a clean non-detection at z — 0.11487, 
but a marginal detection at z — 0.11586. The Li-like ions are a different matter: OVI and NV absorption is detected at both 
redshifts in both lines of both ions (OVI A1038 line at z — 0.11586 is blended with a weak Galactic CI line and is taken as an 
upper limit). Both systems show weak CIV 1548A detections. CIV 1550A is detected for z — 0.11487, but is redshifted out of 
the data range for z = 0.11586. 
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TABLE 2 
Sample HI Absorbers 



Sight Line 


-^abs 


Wr(Lya) 
(mA) 


W,(Ly/3) 
(mA) 


W,(Ly7) 
(mA) 


W,(Ly<5) 
(mA) 


bm 

(km s"-"-) 


log A^Hi 
{N in cm~^) 


method^ 


Mrk335 
Mrk335 
Mrk335 


0.00653 
0.00759 
0.02086 


210 ± 19 
78 ± 3 
149 ±4 


59 ±8 
< 13 
42 ± 7 


< 205 


< 29 


20tl 
33 ±4 
60 ±2 


13.99+0^; 

in 24+0.06 

13.51 ±0.02 


COG 
Ly« AOD 
Lya AOD 



HE 0226-4410 
HE 0226-4410 
HE 0226-4410 
HE 0226-4410 
HE 0226-4410 
HE 0226-4410 
HE 0226-4410 



0.01216 
0.02682 
0.02729 
0.04129 
0.04609 
0.05177 
0.05814 



50 ± 3 
66 ± 10 

30 ± 5 

40 ±4 
157 ± 12 

26 ±3 
30± 10 



< 18 
23 ±4 

< 18 

< 19 
41 ±3 

< 26 



< 14 



< 29 



41 ± 7 

30 ±4 

31 ± 5 

32 ± 7 



16 



+4 
-3 

17 



13 06+" °* 

1 2 7q+0.13 

12.93tH? 
13.841^^9 

1 2 7K+0-27 



Ly« AOD 
Lya AOD 
Ly« AOD 
Lya AOD 

COG 
Lya AOD 
12 ±3 12.77 ±0.14 Lyg fit 



"HI solution method: AOD=integrated apparent optical depth, 
solution, COGt=poorly-constrained curve of growth. 



fit=Voigt fit, COG=multi-line curve of growth 



TABLE 3 



IGM Diagnostic Lines 


ion 


Arost (A) 




[X/H]q^ 


I.P.'^(eV) 


/ion,CIE 


log Tniax 


HI 


1215.67 


0.4164 


0.00 


<13.6 




<4.19'* 




1025.72 


0.07914 










Felll 


1122.52 


0.0544 


-4.55 


16-31 


0.89 


4.27-4.58<^ 


Si HI 


1206.50 


1.63 


-4.49 


16-35 


0.90 


4.30-4.85<i 


CHI 


977.02 


0.757 


-3.61 


24-48 


0.83 


4.66-5.01<i 


Si IV 


1393.76 


0.514 


-4.49 


35-48 


0.35 


4.8° 




1402.70 


0.255 










CIV 


1548.20 


0.191 


-3.61 


48-65 


0.29 


5.0° 




1550.77 


0.095 










NV 


1238.82 


0.157 


-4.22 


78-98 


0.24 


5.25° 




1242.80 


0.078 










OVI 


1031.93 


0.1325 


-3.34 


114-138 


0.22 


5.45° 




1037.62 


0.0658 











Line oscillator strength (Morton 2003)*^ Solar photospheric element abun- 
dance log(X/H); [Xs plund et al. (2005)° Peak ion fraction under CIE via 
[Sutherland fc Dopital ljl993l '): value used in flion and Zion calculations over the listed 
range in ionization potential.'* Tempera ture range over which ion is dominant un- 
der coUisional ionization equilibrium, via lSutherland &: Dopital l|1993l )° Ion is never 
dominant under CIE. Listed temperature is peak CIE abundance. 
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0.29 at logTniax = 5.00. Non-equilibrium ionization 
effects may alter th e post-shock ionization fractions 
(jRaian &: Shulll [20081 ). Because N V lies intermediate in 
ionization potential between the two other Li- like ions, it 
should be a good WHIM tracer as well. However, the so- 
lar nitrogen abundance is more than five times lower than 
oxygen, and nitrogen is often observed to be underabun- 
dant in Galactic HVCs (jGibson et al.ll200H[Coriins et all 
[200l . Until now, N V absorption is largely an unknown 
quantity in the low-z IGM. 

The UV lines of Si IV AA1393.755, 1402.770, Si III 
A1206.500, C III A 977.020, and Fe III A1122.526 are aU 
strong transitions of abundant metals, with production 
ionization potentials below 3 ryd (40.8 eV) and expected 
to probe photoionized, metal-enriched material. In these 
lines, we also have two sets of adjacent ionization levels 
(CIII/IV and Silll/IV), which may help to define the 
ionization state of the absorbers, independent of elemen- 
tal abundance ratios. Highly ionized metals were the pri- 
mary focus of this survey and singly-ionized species were 
not explicitly studied. Howeve r, low - ion absorption ex- 
ists in some Lya IG M systems ([Tripp et al. 2008) and in 
many ionized HVCs ([Sembach et al.lll99& : iCollins et al.l 
120031 '2007). In a future paper, we will examine the be- 
havior of Si II, C II, Fe II, and similar species in the IGM 
and give further analysis options in adjacent ionization 
stages (Sill/III/IV, CII/III/IV, FeH/IH). 

We examined all 13 transitions for each IGM absorber; 
Figures 1 and 2 show examples. In each case, we mea- 
sured the line using the apparent optical depth (ADD) 
method (Sembach & Savage 1992) and/or fitted the line 
with a Voigt profile (see discussion in Paper 2 and sources 
therein). While the Voigt and ADD methods generate 
equivalent line measurements and uncertainties in most 
cases, each has its advantages; ADD is preferable for 
weak, asymmetric, or noisy lines, while Voigt fits are 
more reliable for saturated or blended lines. In cases 
where no line was apparent, we measured a 4cr upper 
limit on equivalent width, Wx, according to 

4 A 

" R iS/N)x ' 

where we adopted instrumental resolutions, R — A/AA, 
for the two spectrographs of i?FUSE = 15, 000 and 
^STis/Ei40M = 42,000. In cases where the measured 
W\ < Wmin, we declared the line to be an upper hmit. 
Heavily blended lines, in which no information was avail- 
able, and lines that were redshifted out of the data range 
were noted but ignored. Measurements of less severely 
blended lines and lines with ambiguous identification 
were measured and taken to be upper limits. 

In the second stage of processing, we checked our de- 
tections for consistency. We required the Ly/3 detections 
or upper limits to be consistent with Lya measurements, 
and the metal doublets (C IV, N V, O VI, Si IV) to be 
consistent within the error bars in logiV. The singlet 
lines (C III, Si III, Fe III) were harder to verify this way, 
but we were able to confirm or rule out many based on 
detections or non-detections in similar ions in the same 
absorption system. 

We carried out a third-stage vetting process of all 
metal-line de tections to elim inate likely fixed pattern de- 
tector noise ([Sahnow et al.|[2000 ) masquerading as weak 
detections, instrumental artifacts, and other inconsisten- 



cies. In some cases a "detection" was indistinguishable 
from adjacent data fluctuations and there was no other 
supporting evidence such as detections in other ions in 
the same absorber. Such cases were deemed suspect and 
demoted to upper limits even though they were techni- 
cally > 4(7. The effects of this process on the detection 
statistics varied. Unfortunately, an unknown number of 
weak but bona-fide detections may have been eliminated 
in this check, but the remaining detections are much 
more robust than in the earlier stages. The elimination 
of some marginal detections is apparent in the low-A^ 
turn-overs seen in some dM / dz plots (see below for dis- 
cussion) . If we assume that absorbers follow a power-law 
distribution with iV, we can assess the effect of this final 
detection verification step. However, the low-A^ distribu- 
tion of metal-line absorbers may also be truncated owing 
to the l imited volumes over whic h galaxies can disperse 
metals ([Tumlinson fc Fan3l2005l Paper 1). 

The final processing stage involved deriving column 
densities and 6-values for the different species. We de- 
termined most (69%) of the H I values directly from Lya 
measurements. Paper 2 showed that Lya-only measure- 
ments underestimate A^hi, owing to saturation and the 
lack of CoG information from higher Lyman lines. The 
effect is strongest in saturated lines and is probably neg- 
ligible for Wlyq ^5 100 tclK. The median equivalent width 
in our sample is W^ya = 94 mA. Eleven H I absorbers 
were measured from Ly/3 only, when Lya was strongly 
blended or redshifted out of the E140M spectral range 
{z > 0.4). Values of column density and doppler width 
for strong H I absorbers were usually (30%) determined 
from concordance CoG plots (see Paper 2 for discussion) . 
In some cases, we measured Ly7 and hyS to further con- 
strain the H I solution. Lyman-line equivalent widths 
and the derived 6hi and A^hi measurements are listed in 
Table [2 

We took column densities and 6-values of the doublet 
ions (O VI, N V, C IV, Si IV) as the error-weighted mean 
of the two lines, when possible, or as measured if only 
one line was available. Values for the single lines (C III, 
Si III, Fe III) were taken as measured by the AOD or 
Voigt profile methods. We present truncated tables of 
the measured equivalent widths, 5-values, and column 
densities for all seven ions in Tables 4-10. Note that 
AOD-measured 6-values are equivalent to Voigt fit mea- 
surements for strong lines, but become increasingly un- 
certain for weak and noisy lines. Some of the poorly- 
determined 6-values are listed as approximate and AOD 
5-values in the tables should generally be treated with 
caution. Full versions of the Tables, including both de- 
tections and upper limits, are available in the electronic 
edition. 

2.4. Absorber Statistics 

With such a large database of detections, we are able 
to analyze the absorber statistics with considerably more 
confidence. We constructed histograms showing the 
number of absorbers, A/i, in logarithmic bins of width 
0.2 dex in column density. In general, the median column 
density uncertainty. Slog N, for the metal ion detections 
is larger than half a bin size. We analyzed the distri- 
bution A/'(log A^) using a "distributed-power" histogram 
technique to spread the power of each detection across 
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TABLE 4 
Sample O VI Detections 



oignt Line 




\LiJo2) 




L / 1 AQ Q^ 

(lOooJ 




0[U Vlj 


lOgiV^U V L) 


notes 






(km s ^) 


(mA) 


/I — 1 ^ 
(km s ^) 


(mA) 


/I — 1 \ 
(km s ^) 


(N m cm *') 




HE 0226-4410 


0.20700 


34 ± 1 


184 ± 10 


34 ±8 


111 ±4 


34 ±8 


14.33 ±0.06 


fit. AOD 


HE 0226-4410 


0.22004 


~ 16 


12 ± 1 




< 31 


~ 16 


13.05 ±0.15 


AOD 


HE 0226-4410 


0.24515 


~ 29 


23 ± 1 


~ 18 


10 ± 1 


~ 25 


13.29 ±0.15 


AOD, AOD 


HE 0226-4410 


0.31833 


6±3 


9±2 




< 7 


6 ± 3 


12.97 ±0.17 


AOD 


HE 0226-4410 


0.34031 


21 ± 2 


74 ±11 




< 82 


21 ± 2 


13.75 ±0.05 


fit 


HE 0226-4410 


0.35493 


13 ±9 


9±3 




< 8 


13 ±9 


12.94 ±0.26 


AOD 


HE 0226-4410 


0.35524 


44 ±6 


63 ± 2 


35 ± 5 


22 ± 1 


41 ±6 


13.69 ±0.08 


AOD, AOD 


PKS 0312-770 


0.20192 






~ 40 


31 ± 1 


~ 40 


13.76 ±0.15 


AOD 


PKS 0312-770 


0.20275 


62 ± 2 


527 ± 28 


66 ± 2 


323 ± 16 


63 ±2 


14.85 ±0.07 


fit, fit 


PKS 0405-123 


0.09657 


40 ±20 


59 ± 31 


14 ± 13 


32 ± 27 


31 ±20 


13.71 ±0.21 


fit, fit 


PKS 0405-123 


0.16714 


82 ±4 


470 ± 40 


61 ± 5 


198 ± 18 


75 ± 10 


14.68 ± 0.05 


fit, AOD 



TABLE 5 
Sample C IV Detections 



Sight Line z^bs 6(1548) Wr(1548) 6(1550) Wr(1550) 6(CIV) logAf(CIV) notes 







(km s-i) 


(mA) 


(km s-i) 


(mA) 


(km s-l) 


(AT in cm~2) 




PKS 0312-770 


0.10877 


22 ± 1 


78 ± 11 




< 55 


22 ± 1 


13.52 ±0.25 


AOD 


HS 0624+6907 


0.06348 


25 ± 1 


113 + 3 


~ 27 


44 ±3 


25 ± 1 


13.49 ± 0.09 


AOD, AOD 


HS 0624±6907 


0.07573 


10 + 1 


45 + 4 


12 ± 3 


21 ±6 


10 ±3 


13.06 ±0.07 


fit, fit 


PG0953±414 


0.04960 


17 + 5 


33 + 6 




< 22 


17 ±5 


12.98 ±0.16 


AOD 


PG 0953+414 


0.06808 


11 ± 1 


135 ± 24 


10 ± 1 


93 ± 18 


10 ± 1 


13.83 ±0.09 


fit, fit 


PG 1211+143 


0.05104 


46 ± 1 


260 ±4 


33 ± 2 


143 ±2 


41 ±6 


14.04 ±0.04 


AOD, AOD 


PG 1211+143 


0.06447 


46 ± 5 


73 ± 12 


39 ± 6 


49 + 10 


43 ±6 


13.32 ±0.07 


fit, fit 


PG 1211+143 


0.06493 


24 ±5 


31 ± 5 


~ 46 


38 + 5 


24 + 5 


13.10 ±0.21 


fit, AOD 


Mrk205 


0.00443 


50± 10 


362 ± 14 


55 ±2 


235 + 9 


51 + 10 


14.05 ±0.18 


AOD, AOD 


Mrk 205 


0.02746 




< 131 


~ 17 


53 ± 1 


~ 17 


13.55 ±0.10 


AOD 


Q 1230+0115 


0.04886 


21 ±9 


39 ±1 




< 21 


21 + 9 


13.00 ±0.16 


fit 



TABLE 6 
Sample N V Detections 



Sight Line 




6 (1238) 


Wr(1238) 


6 (1242) 


Wr(1242) 


6(NV) 


logAf(NV) 


notes 






(km s-i) 


(mA) 


(km s-l) 


(mA) 


(km s-i) 


{N in cm~^) 




PKS 0312-770 


0.20192 


8 + 1 


25 + 4 


9±2 


17 ±3 


8±2 


13.12 ±0.07 


fit, fit 


PKS 0312-770 


0.20275 


67 + 4 


133 + 10 


56 ±6 


63 ±7 


63 + 6 


13.74 ±0.04 


fit, fit 


PKS 0405-123 


0.16714 


51 + 5 


143 + 16 




< 155 


51 + 5 


13.83 ±0.04 


fit 


PKS 0405-123 


0.33400 






~ 12 


28 ±1 


~ 12 


13.52 ±0.12 


AOD 


PKS 0405-123 


0.36079 


24 + 8 


46 ± 11 






24 + 8 


13.47 ±0.12 


AOD 


HS 0624+6907 


0.29530 


18 + 5 


30 + 2 


7±3 


11 ± 1 


14 + 5 


13.17 ±0.13 


AOD, AOD 


PG 0953+414 


0.06808 


8+1 


45 ±4 


15 ± 2 


30 ±5 


10 + 3 


13.46 ±0.05 


fit, fit 


PG 0953+414 


0.23350 


9 + 3 


11 ± 5 




< 8 


9 + 3 


12.68 ±0.14 


fit 


3C 249.1 


0.11925 


~ 15 


11 ± 2 




< 10 


~ 15 


12.78 ±0.20 


AOD 


3C 249.1 


0.17371 


24 + 7 


29 ±4 




< 17 


24 + 7 


13.19 ±0.13 


AOD 


PG 1116+215 


0.05904 


~ 27 


20 ±2 


~ 19 


18 ±2 


~ 24 


13.16 ±0.12 


AOD, AOD 



several log Ni bins. The total distribution is then, 
Af{logN) = Y,mogNi) (2) 

i 

sr Q-^ /-[iog7v_-bg7v^\ 

4" Slog N,V2^ '''''' I 2 [Slog iV,] 2 )■ 

The difference between this distributcd-power method 
and simple histogram binning is subtle; sharp peaks and 
valleys in the simple histogram are smoothed to some 
degree, but the general shape remains unchanged. We 
believe this method takes measurement uncertainty into 
better account. 
At low column densities we must correct the survey 



for incompleteness (Papers 1 and 2). To do so, we divide 

the number of absorbers of a given equivalent width W 
by an "effective redshift pathlength" , Azcs(W), which 
accounts for absorption-line sensitivity. For each species, 
we calculate the number of absorbers per unit redshift, 
dJ\f /dz, as a function of column density, N, to account 
for incompleteness in the weaker absorbers, 

dAA(log N) _^ mogNi) 
I 

We measured the sensitivity by the signal-to-noise ra- 
tio, S/N, determined for each sight line as the ratio of 
mean flux to standard deviation in the flux in the line- 
free continuum of each 10 A segment of FUSE or STIS 
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TABLE 7 
Sample C III Detections 



bight Line 




Wr(y/ I) 


V.t TTT\ 
0(U 111 ) 


log JV [K^ llij 


notes 






(mA) 


/I — 1 \ 
(km s ) 


(N in cm 




HE 0226-4410 


0.20700 


219 ± 66 


23 + 5 


13.98 + 0.24 


fit 


HE 0226-4410 


0.34031 


23 ±3 


32 + 4 


12.63 + 0.10 


ADD 


PKS 0312-770 


0.20192 


270 ± 17 


57 + 4 


13.63 + 0.20 


ADD 


PKS 0312-770 


0.20275 


600 ± 35 


94 + 5 


13.97 + 0.35 


AOD 


PKS 0405-123 


0.16714 


369 ± 73 


16 + 2 


14.78 + 0.71 


fit 


PKS 0405-123 


0.33400 


36 ±4 


38 + 2 


12.82 + 0.07 


AOD 


PKS 0405-123 


0.36079 


141 ± 10 


17 + 1 


13.61 + 0.05 


fit 


HS 0624+6907 


0.32088 


10 ±3 


12 + 4 


12.10 + 0.13 


fit 


PG 0953+414 


0.06808 


112 + 20 


28 + 4 


13.35 + 0.06 


fit 


PG 0953+414 


0.11827 


36 + 3 


26 + 8 


12.83 + 0.09 


AOD 


PG 0953+414 


0.14232 


30 + 16 


20 + 12 


12.65 + 0.19 


fit 



TABLE 8 
Sample Si IV Detections 



Sight Line 




b (1394) 
(km s"-') 


Wr(1394) 
(mA) 


b (1403) 
(km s~^) 


Wr(1403) 

(mA) 


b(Si IV) 
(km s~^) 


logAf(Si IV) 
(A^ in cm~^) 


notes 


PKS 0312-770 


0.20192 


35 + 8 


177 + 54 


40 + 2 


80 + 7 


36 + 8 


13.39 + 0.11 


fit, AOD 


PKS 0312-770 


0.20275 


76 + 3 


595 + 99 


75 + 1 


470 + 12 


75 + 3 


14.17 + 0.19 


fit, AOD 


PKS 0405-123 


0.16714 


33 + 2 


164 + 12 


25 + 3 


117 + 21 


30 + 4 


13.38 + 0.09 


fit, fit 


HS 0624+6907 


0.06348 


24 + 2 


47 + 4 




< 54 


24 + 2 


12.80 + 0.10 


AOD 


PG 1116+215 


0.13850 


17+ 1 


34 + 2 




< 30 


17 + 1 


12.69 + 0.08 


AOD 


PG 1211+143 


0.05104 


15 + 1 


73 + 5 


17+ 1 


45 + 4 


15 + 1 


13.02 + 0.03 


fit, fit 


Mrk205 


0.00443 


30+ 1 


115 + 9 


24 + 9 


74 + 7 


28 + 9 


13.30 + 0.06 


AOD, AOD 


Q 1230+0115 


0.07805 


~ 10 


29 + 4 


~ 20 


24 + 4 


~ 13 


12.70 + 0.16 


AOD, AOD 


PG 1259+593 


0.23951 


10 + 8 


37 + 30 






10 + 8 


12.61 + 0.29 


fit 


PKS 1302-102 


0.04225 


~ 18 


18 + 2 




< 45 


~ 18 


12.39 + 0.23 


AOD 


PKS 1302-102 


0.09485 


~ 18 


18 + 7 




< 12 


~ 18 


12.36 + 0.27 


AOD 



TABLE 9 
Sample Si III Detections 



Sight Line 


-^abs 


Wr(1207) 


6(Sini) 


logAr(Silll) 


notes 






(mA) 


(km s'l) 


(A'' in cm~^) 




HE 0226-4410 


0.16340 


19 + 1 


~ 9 


12.07 + 0.09 


AOD 


HE 0226-4410 


0.19859 


20 + 3 


26 + 3 


12.07 + 0.12 


AOD 


HE 0226-4410 


0.20700 


49 + 3 


24 + 4 


12.47 + 0.05 


AOD 


HE 0226-4410 


0.30565 


30 + 7 


12 + 7 


12.29 + 0.11 


AOD 


PKS 0312-770 


0.05946 


21 + 5 


21 + 6 


12.07 + 0.15 


AOD 


PKS 0312-770 


0.17976 


13 + 8 


8 + 5 


11.77 + 0.21 


fit 


PKS 0312-770 


0.20192 


123 + 5 


30 + 1 


12.97 + 0.05 


AOD 


PKS 0312-770 


0.20275 


647+ 151 


36 + 2 


14.92 + 0.27 


fit 


PKS 0405-123 


0.16714 


266 + 12 


30 + 1 


13.34 + 0.02 


fit 


PKS 0405-123 


0.21070 


104 + 12 


40 + 9 


12.86 + 0.05 


AOD 


PKS 0405-123 


0.35210 


128 + 17 


36 + 4 


12.74 + 0.05 


fit 



data. In regions of overlapping data (e.g., the overlap 
of FUSE LiFl and LiF2 channels or the overlap between 
FUSE and STIS/E140M coverage), we used the higher 
S/N value. The vector, {S/N)\, for a given wavelength 
A was rebinned to a resolution comparable to that of the 
data. 

We then modified {S/N)\ to take into account regions 
in which IGM absorption lines would not be detected. 
We used two methods to assess the sensitivity of the ef- 
fective pathlength, Az, to different {S/N)\ algorithms. 
In the first method, we set {S/N)\ = in ~ 50 re- 
gions of strong interstellar lines and instrumental arti- 
facts that appear in most or all of the data sets. In the 
second method, we let the smoothed, normalized flux 



data mask the {S/N)\ profile. Regions of strong absorp- 
tion {t\ ^ 1) were completely masked, {S/N)\ — 0, 
while areas of unsaturated absorption were masked pro- 
portional to the optical depth, {S/N)x oc (r^ -I- 
This effectively removed fixed ISM lines and instrumen- 
tal artifacts, as well as strong IGM lines and intrinsic 
AGN features. It also accounted for degraded sensitivity 
in regions with moderately strong absorption lines of any 
type (AGN, IGM, or ISM). 

Next, we converted each {S/N)\ vector to a 4ct mini- 
mum equivalent-width vector using equation [T] For sin- 
gle lines, we determined the total redshift pathlength, 
Az{W), at a given equivalent width sensitivity by sum- 
ming the available pathlength over which W\ > Wmin 
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TABLE 10 
Sample FeIII Detections 



Sight Line 


^abs 




J,/T?^ TTT^ 

o(re ilij 


log Jv (re 111 ) 


notes 








(mA) 


(kra s ) 


/ JV7- ■ O \ 

(N m cm ) 




HE 0226- 


4410 


0.16236 


23 ±4 


~ 20 


13.64 ±0.14 


AOD 


PKS0312 


-770 


0.20192 


20 ±7 


6±2 


13.52 ± 0.11 


fit 


PKS0312 


-770 


0.20275 


127 ±5 


64 ± 1 


14.46 ± 0.04 


AOD 


PKS 0405 


-123 


0.16714 


22 ±5 


10 ±2 


13.56 ± 0.09 


fit 


PKS 0405 


-123 


0.35210 


28± 15 


16 ±6 


13.58 ± 0.17 


fit 


3C 249.1 




0.10834 


15 ±1 


~ 17 


13.47 ±0.22 


AOD 


3C 249.1 




0.14183 


36 ±4 


19 ±5 


13.85 ± 0.08 


AOD 


3C 249.1 




0.24396 


16 ±3 


~ 4 


13.63 ± 0.12 


AOD 


PG1116+ 


215 


0.13850 


14 ±1 


~ 9 


13.44 ±0.10 


AOD 


PKS 1302 


-102 


0.07695 


24 ±3 


~ 18 


13.69 ± 0.13 


AOD 


Mrk876 




0.11487 


25 ±7 


16 ± 5 


13.60 ±0.11 


fit 



in all 28 sight lines. For multiplet lines, including Lya 
and Ly/3, we scaled a second vector, Wmin{z), by the 
line-strength ratio, /A, and shifted to the rest frame of 
the weaker line. In this case, Az{W) is the sum of all 
available pathlength for which either W\ > Winin(2i) 
or Wx > [Wmin(2:2)](/iAi//2A2). In all cases, we ex- 
cluded pathlength within Az < 0.005 (1500 km s"^) 
of the intrinsic AGN redshift and within Az < 0.0017 
(500 km s-^) of the Local Standard of Rest (LSR). Thus, 
systems intrinsic to the AGN (outflows) or to our Galaxy 
( high velocity clouds ) were not included. 

iTripp et al.l (2008) study both intrinsic and interven- 
ing systems and adopt a more restrictive velocity limit 
(Az > 0.017, V > 5000 km s^^) than we use here. 
As they note, there are numerous examples of high-z 
quasars where intrinsic absorption components are seen 
at V > 5000 km s~^. They find the line density for intrin- 
sic absorbers is roughly a factor of two higher than for 
intervening systems. Nevertheless, Tripp et al. note that 
the excess of systems is largely confined to the velocity 
range —1000 < t'dispi ^ 2500 km and that intrinsic 
systems are relatively rare at Wdispi > 2500 km a~^. 

Both methods of {S/N)\ masking reduced the total 
Az in the data by 3-4% relative to the unmasked S/N 
profile. However, the difference in total Az between these 
two methods was less than 1%. The second, self-masking 
method allocates relatively less pathlength for weak lines 
and relatively more for stronger lines. We deemed this 
method to be more realistic and used it throughout our 
analysis. 

In our previous work, the absorber redshifts were gen- 
erally small and the redshift is linearly proportional to 
comoving distance, di = (c/Hq) dz. In the new survey, 
we included absorbers out to redshifts z = 0.4, where a 
given a linear distance d£ covers a larger redshift interval 
dz than at z = 0. We therefore use the relations. 



z)"^ dz, so that 



di _ 1 c 
d^^ (l + z) H{z) 
c 1 



Ho (1 



z 



[n^{i + zf + n^\-^/^ 



(4) 



where we assumed a flat (rim, A) cosmology with H{z) — 
Ho [n^{l + zf + n^]^/'^. By including the (1 + z)^ den- 
sity scaling and normal izing to Hq, it is conventional 
(jBahcall fc Peeblei Il969i) to define the dimensionless 
"absorption pathlength function", dX = [Ho/H{z)]{l + 



dX={l + zf [ilrail + zf + ilAr^^'^ dz 



(5) 



We calculated AX{Ni^in) from the individual spectra 
in the same manner as Az^N^i^) above using numerical 
integration and ad opting 0^ = (0-261 ± 0.016)hjQ and 



f^A = 0.716±0.055 (ISpereel et al.ll200l . For low-redshift 
absorbers, such as all C IV detections (zabs < 0.11), the 
difference between Az and AX is negligible. We list 
both dN/dz and dN/dX in the discussion, to provide 
both observational and theoretical perspectives. 

Uncertainties in dM /dz arise as uncertainties in both 
N a nd Az. In the numerator, one-sided Poisson statis- 
tics (lGehrelslll986h for dNi are assumed. While the un- 
certainty in Az(iy) is negligible as discussed above,, we 
must still take into account the changing Az over the 
width of a column density bin. For weak lines where the 
Azoff(M^) is falling sharply, there is a significant differ- 
ence in Az between one side of an equivalent width or 
column density bin and the other. We approximate this 
uncertainty as dAzi — i[(Az)i+i — (Az)i_i] or roughly 
the change in Az over a quarter of a bin, the typi- 
cal Icr dispersion of column densities from bin center. 
The dM I dz uncertainties from numerator and denomi- 
nator are added in quadrature. Uncertainty from Pois- 
son statistics dominates in almost all cases, but dAzi 
becomes significant for numerous absorbers (H I, O VI, 
Si III) in the weakest column density bins. Uncertainties 
in dM jdX are completely analogous. 

2.5. Baryon Fraction 

Through measurements of D /H and acoustic peaks in 
the CMB, we now have an accurate estimate of the cos- 
mological density o f baryons, ^70 = 0.0455 ± 0.0015 
(jSpergel et al.|[2007l ). However, most of these baryons do 
not exist in collapsed form, and probabl y are distributed 
throu ghout various phases of the IGM (jCen fc Ostrikeil 
119991: IShullll2003f ). This is certainly the case at high red- 
shift, where the Lya forest is the repository for large 
quantities of gas that has not yet collapsed into galax- 
ies. But even at low redshifts, the IGM likely contains 
between 50% and 80% of the baryonic matter. Because 
low-density gas at a variety of temperatures can be dif- 
ficult to detect, a complete cosmological baryon census 
requires a combination of ultraviolet and X-ray spectro- 
graphs. 

Here, we use UV spectroscopic column densities and 
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absorber frequencies of H I and various metal ions to de- 
rive the baryon density contained in IGM thermal phases 
at T < 10^ K. To compare to similar estimates in the lit- 
erature, it is convenient to define two such measures. The 
first, directly measured quantity provides the contribu- 
tion to the critical mass density, of the IGM absorber 
mass in a given ion, 

f^ion = I I X (6) 



X 

W(log Mon) 



N„ 



dX 



The integration is over the column density distribution 
and the corrected redshift pathlength parameter dX. 
Here, w-ion is the ion mass, and the critical density is 
p„ o = (3iJ^/87rG) = (1.879 x IQ-^s g cm~3)/i2^ The 
integration is over the column density distribution and 
the corrected redshift pathlength parameter dX. 

Uncertainty in fiion is taken as the asymmetric Pois- 
son uncertainty on dN{> N) evaluated at a d/V- 
weighted mean dX, (z), and {N). The value of fJion oc 
[Z (M/7J)0 /ion/i7o]^^7 so we do not include systematic 
uncertainties in these quantities. Non-equilibrium ion- 
ization may produce higher ion fractions for Si IV and 
the Li-like ions while photoionization may lower the as- 
sumed near-unity CIE ion fractions of C III, Si III, and 
Fe III. Similarly non-solar abundance ratios or deviations 
from the canonical 10% solar metallicity often assumed 
for IGM clouds will introduce systematic differences. 

The second parameter, fijcM' estimates the density of 
all baryons traced by the particular ion after correcting 
for abundances and ionization stages. In Paper 1, we cal- 
culated the fraction of the universe composed of baryons 
at T = 10^~^ K, using O VI as a tracer for baryonic 
WHIM gas in general. This parameter includes correc- 
tions for the metallicity, Z , of the element relative to 
solar abundances, and for the ionization fraction, fovi, 
of the particular ion. Our adopted formula was 

(OVi) _ ( Hp \ firUH ^ 



a 



IGM 



cpci fij Z {0/H)q fovi 
'dAf{log NoYi) 



N„ 



dz 



{Novi) d\og NoYi , 



where fi — 1.32 for low-meta llicity gas and primo rdial 
helium, Yp « 0.2477 ± 0.0029 (jPeimbert et al.ll2007D . 

We can easily generalize this formula for any ion divid- 
ing flion by an assumed metallicity, Z = 0.1 Zq, abun- 
dance of element M relative to hydrogen, {M/H)q, and 
ionization fraction (assumed to be the peak CIE ion 
abundance /ion listed in Table [3]). Because the redshift 
range in our current work is larger than that covered in 
Paper 1, we adopt pathlength coordinate X instead of z, 
as discussed above, and we perform the sum over bins of 
equal AX. 
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ion / 1 J 



A log A^ionAXj- 



x„ 



E E 

j = i = -/Vmi: 

As in our previous work, we adopt a fiducial IGM metal- 
licity of Z = 0.1^0 and a Hubble constant Hq = 
(70km s"^ Mpc"^)/i7o. We sum dM/dX for each ion 
over bins in column density Mon and normalized path- 
length dX . The final term, AX^/Xmax, averages the 
dJ\f/dX contributions from each dX bin. 

2.5.1. Baryons in the Lya Forest 

While the main focus of this paper is metal ions, our 
database of H I absorbers is the largest to date in the 
low-z IGM and can provide good constraints on S^jq^, 
the fraction of IGM bar yons in the Lya forest. For this 
we use the metho ds of jPenton. Shull. fc Stocki (|2000l ) 
and ISchavi ()2001[ ). which are slightly different from the 

methods employed for rij^J^j^ above. In the first method, 
the IGM absorbers are assumed to be isothermal spheres, 
and the total cloud mass is inferred from the H I column 
density and an assumed impact parameter. The neutral 
fractio n is derived from photoionization conditions. Fol- 
lowing |Pe^on;iSh^iL&_Stockc (2000 ), we define 
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x=a log N^i 



dA/'(log 



dz 



Nm) ^1/2 



AX 



"14 



(9) 



where ^loo is the typical scale of Lya forest clouds 
scaled to 100 kpc, J-23 is the ionizing ra diation field 
intensity at low redshift (jShull et al.l Il999f ) in units of 
10"^'^ ergs cm~^ s~^ Hz~^ sr"-*^, a^ is the spectral index 
of the radiation field above 1 ryd scaled to as = 1.8, 
and iVi4 is the H I column density in units of 10^ "^ cm"^. 
Secondly, we use the formalism of ISchavj (|2001h . which 
assumes that IGM clouds are gravitationally bound and 
that their observed column densities are typical over a 
Jeans length: 



(HI) 
IGM 



(1.46 X io-4)/i7o'rL^i'2?;°-'''x 



N„ 



4«ll = (1.46xlO-^);^,-l^L/f,Tr^x 

^max log iV„a, 



E E 

X=0 log W„,: 



dAf{log TV] 
dz 



mj ^1/3 



AX 



(10) 



As with fiionj we perform a two-dimensional sum 
over both column density and redshift bins assuming 
r = 2 X lO'* K in the photoionizcd gas and a hy- 
drogen photionization rate Fhi = 0.03 x 10^^^ 
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assumptions about ionization corrections, to infer a sta- 

(see the derivations and 
We quote each me tallic- 



TABLE 11 tistical metalhcity of the IGM (see the derivations and 
IGM Detections and ResuliDS Summary ti ' i j o^ iir i u j. n- 
discussion 111 Papers 1 and zj. — — 



^abs 



ai4 



Cl4 



. \t>y «is a tWpUOia o: 

HurfWo thcni 



ar abundance J Aspiund et al.l 
y Jioii: the peak ion abunda nce m 



OVI 

NV 
CIV 

cm 

Si IV 
Si III 
Felll 
HI 



83 
24 
24 
39 
20 
53 
14 
650 



< 0.40 

< 0.396 

< 0.116 

< 0.40 

< 0.24 

< 0.40 

< 0.40 

< 0.40 



5.22 
5.30 
2.42 
4.84 
4.21 
5.14 
5.11 
5.27 



0.71 ± 0.03 
0.98 ± 0.07 
0.87 ±0.06 
0.56 ± 0.05 
0.74 ±0.06 
0.81 ± 0.03 
0.71 ± 0.05 



2.0 ±0.1 
5.3 ±0.7 
5.5 ±0.7 

9.5 ±0.9 
12 ±2 
41 ±3 

1.6 ±0.2 



O.l5cbni8ional i5i|zationlfei|iilibri.i9ftti(SUtherlalid & Dopital 
0.3£']dg9'§Ji or ia^hotoiofiite|tion%^MiWaim (sfee Tabled. 
0.055^^0g5fgj.j.(lg_22etaUi(b?ti|s arfi-1&^^fe?age ffer the entire 
0.0l5g^^^^. -vir^q^jjg no'^S^inis metalhcity of in- 

0.l9(|ljT^g^al absorbers. 'Bm obserV-^d^^dWibly-iokized species 
0.023(:^0i$f2si ii5+|^e III) Save peA8€*B-a%und£ftices close to 

^o^nf (C iV, N V, OVI) 



Previous Work/Literatuifea^fluBgak CIE ionization tractions /peak U.2 - 0.3 

^ ^|n coluny], 8. we provide thc^^Cjt/il^ aboorljpr fr equency 
per ""'^ '''' ' 



unit r§^t|^ift, dSfJdz, 
equivalent8^si?|t§i 30 mA (the survey goes4deeper than 
ojijiite limit iiEt|any cases) . Cql.ig8iiio.@4gives yalues of the 
propcr-leiipMfcrrected dW/S^^jOW > <80 mA. The 
_ . , , cumulative bivariate distribution of absorbers per red- 

1 nis Work - equivalent TfaBipIe ^ , ^ ,. r ^ ^ n..i 
shift intprval a.s a tiinrlinn nt pnlnmn ripnsity wa.s titled 

O.l^MBla po^*8r:|aw, [#34^> A2)^gi(H^&K N] k N~^\ and 



OVI 


40 


< 0.15 


2.21 


0.9±0.1 


2.5±0.2 


OVI 


44 


0.12-0.57 


~ 1.9 






OVI 


53 


< 0.5 


3.18 


~ 0.9 


~ 2 


cm 


30 


< 0.21 


2.41 


0.73±0.08 


11.5±1.0'^ 


HI 


187 


< 0.069 


0.770 







infcgrkted dowji to fiducial 



OVI 
CHI 



35 
28 



< 0.15 

< 0.21 



2.96 
3.485 



0.81±0.06 
0.65±0.05 



1.96±0.19 
7.69±0.91 



0.0l6t;fe%q§giex /f^alisted ifttaolumiiOf S:.0.15ie po-^er-law slope 
becomes systematically steeper by . 05-0 . 1 de x when 



becomes systematically steeper by U . Uo-U . i de x when 

^ Scaled by the maximum CIE abundance (Sutherland & Dopita 1^^^^ Makh'Mm\&irmWimiih]Sm%y^^mc of the 

integrated down to equivalent width ly > 30 mA'^ 1-this work; 2-Danfog^|^ S^^lk^'l,?&ii"T8Ifhg* l!igl^'l'&^ the dis- 

2007; 5-Danforth et al. 2006; 6-Penton et al. 2004^^ Published value of pararneter Ci4 (Paper ,,21 and, fiSncfi Zn = 12% wjere 
incorrect. Corrected values are listed here. tflButlOn USm^ 0& distributed pbwer teclinique. We 



(|Shull et all fTOOl iWevmann et~al\ [200lh . Uncer- 
tainties in rijQ|][(log N) are the quadratic sum of 
the one-sided Poisson error for total redshift path- 
length dM{dlogN,dX). iPenton. Stocke. fc ShuU (l200l l 
showed that uncertainty due to cosmic variance is small 
relative to Poisson error for Az > 1. Since we are sum- 
ming over Az > 5 for all but the lowest column density 
bins, we discount cosmic variance as a significant source 
of error. 

3. RESULTS AND DISCUSSION 

In Table llll we present a compilation of our results. 
The first three columns list ion name, total number of 
detections (at any equivalent width), and redshift cov- 
erage of the survey. Column 4 lists the maximum total 
redshift pathlength sampled to any sensitivity; because 
of (S/N) dependence, a smaller total pathlength is sen- 
sitive to weaker absorbers. 

For each ion, we produce a plot of the "multiphase 
ratio", log (TVHi/Mon) vs. logA^ni- This ratio was ini- 
tially defined for O VI (Paper 1) to characterize the po- 
tential connections between different tracers of baryonic 
mass (H I and O VI) observed to be kinematically as- 
sociated in the low-z IGM. We assume that absorbing 
clouds in the IGM have both a photoionized component 
(traced by sharp Lya lines and low-ionization metal ions) 
and shock-heated gas (WHIM, traced by O VI or other 
high-ions). CoUisionally ionized C III, Si III, etc. are 
possible, but we consider this less common due to the 
fast cooling rates at T < 10^ K, where these ions reach 
their peak CIE abundance (see Paper 2 for discussion). 
In columns 5 and 6, we list the slope, ai4, and the in- 
tercept, Ci4, at column density A^hi= 10^'' cm~^ of an 
error-weighted linear fit to the multiphase correlation. 
These two parameters are used, together with various 



determined /3 by an error-weighted fit to the cumulative 
dM /dz. However, the index fitted to cumulative dM /dX 
is comparable, within the uncertainties in all cases. 

Where possible. Table [TT] lists the equivalent values 
from Papers 1 and 2, as well as several other large surveys 
of the low-z IGM. For direct comparison of the current 
survey to our previous work, we calculate O VI statistics 
for z < 0.15 and C III statistics for z < 0.21, consistent 
with the limits in Papers 1 and 2, respectively. We have 
not adopted the Wlyq > 80 mA criterion used in those 
papers, nor have we accounted for the somewhat less 
stringent detection requirements adopted in our earlier 
work (detections in Papers 1 and 2 were 3cr, while the 
current work requires a minimum of 4(t). 

We continue our results in Table [121 listing fiion 
and fijcM for absorbers with equivalent width lim- 
its Wx > 10 mA and Wx > 30 mA. The 17;^'^ 
values are relative to the fiducial b aryon fraction 
^^baryon = 0.0455/iyQ^ (jSpergel et al.ll200"7[ ) and are scaled 
by (/pcak//ion) {Z/0.1)-'^ h^Q as described above. As in 
Table [Til above, we list related values from previous sur- 
veys and compare to equivalent values from previous sur- 
veys. 

3.1. Discussion of Individual Species 
3.1.1. HI 

Our new H I r esults (Figure [3l) are compatib le with 
previous studies. IPenton. StockeT fc ShulU (|2004( l found 



dM/dz 



165 ± 15 for Lya absorbers down to Wx 



30 mA equivalent width, with a median VFlyq — 68 
mA and a power-law distribution in column density 
dAf{N)/dz cx N-'^ with /3hi = 1-65 ± 0.07. In Paper 2, 
we searched for O VI and C III absorbers in Lya systems 
with Whya > 80 mA, correcting for Lya saturation with 
full CoG determinations of column density in most cases. 
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log N(HI) log N(H1) 

Fig. 3. — HI detection statistics and dM /dz. Left panel shows the number of HI absorbers per logarithmic bin, Alog TVjji = 0.2, 
uncorrected for completeness. Note the turnover at logTVjji < 13. Dashed curve shows absorbers for which higher-order Lyman lines were 
used to constrain the (A'^hIi ''hi) solution. Dotted curve shows Lyo-only solutions. Right panel shows c ompleteness-corrected differential 
dj^ /dz vs. A'^HI- Error bars on the left panel are based on one-sided Poisson counting errors IIGehrels|[l986i) . while the error in dLM /dz 
includes a component of Az uncertainty as well. See text for discussion. 



TABLE 12 
Oion AND fijQM Summary 



ion 


f^ion [10-8] 


fiion [10-«] 


''iGM/^'b 


^iGM/f^b^ Source^ 




(> 10 mA) 


(> 30 mA) 


(> 10 mA) 


(> 30 mA) 



OVI 


49 ±4 


41 ± 5 


0.086 ± 0.008 


0.073 ± 0.008 


1 


N V 


3.6 ±0.8 


3.0 ±0.9 


0.050 ±0.011 


0.042 ±0.012 


1 


CIV 


8.2 ± 1.4 


7.7 ± 1.5 


0.027 ±0.005 


0.026 ± 0.005 


1 


cm 


8.9 ± 1.3 


7.8 ± 1.4 


0.010 ±0.002 


0.009 ± 0.002 


1 


Si IV 


8.4 ± 1.1 


8.3 ± 1.1 


0.075 ±0.010 


0.075 ±0.010 


1 


Si III 


11 ± 2 


10 ±2 


0.037 ±0.006 


0.037 ±0.007 


1 


FoIII 


20 ±5 


11 ± 7 


0.040 ±0.011 


0.022 ±0.015 


1 



Previous Work/Literature Values 


OVI {z 


< 0.15) 


26 ± 4 


0.048±0.007 


2 


OVI {z 


> 0.12) 




0.059 


3 


This Work - Equivalent Sample 


OVI (z 


< 0.15) 


49 ±4 33 ±5 0.087 ±0.008 


0.058 ± 0.008 


1 



^ Scaled by fiant the maximum CIE abundance (Sutherland &; Dopita 1993), Z = 0.1 Zq, 
and Ho = 70 km s"! Mpc"!. = 0.0455 ± 0.0015^^ l lSpergel et al.|[2007^ . See 
text.'' 1-this work; 2-Danforth & ShuU 2005; 3- Tripp et al. 2005 



Despite this correction for artificially low values of A^hi, 
we found only a minor change in index, /3hi = 1.68±0.11. 
In our new work, we do not use a minimum Lya equiva- 
lent width, but rather are limited only by the quality of 
the data. We find a survey median of Whya = 94 mA, 
with the uncorrected distribution peaked at Whya ^ 
45 mA. The power-law index is /3hi = 1-73 ± 0.04, 
and the distribution has an integrated line frequency of 
dU /dz = 129tt down to 30 mA, significantly different 
from the result of lPenton. Stocke. fc ShiM (|2004[ ). How- 
ever, all the absorbers in the Penton et al. surveys were 
at z < 0.07, while those in the current survey are at 



z < 0.4 (with (z) ~ 0.14). The proper-length-corrected 
absorber frequency is d/V /dX = 151^g, more in line with 
the lPenton. Stocke. fc Sluilil (|2004[ ) resuh. 

The A'hi histogram and d/V / dz plot show a small dip at 
log A'hi — 13.6 in an otherwise smooth and well-behaved 
distribution. This is likely a selection effect of how we 
chose and measured Lya absorbers. In Paper 2, we 
showed that H I measurements using only the Lya line 
tend to underestimate A'hi because of line saturation. 
The dip at logA^ni = 13.6 lies at the transition between 
Lya-only and CoG-determined column densities and is 
almost certainly an artifact. Since the total number of 
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Fig. 4. — Distribution of doppler b- values for 650 HI absorbers. 
Dashed line represents the subset of 195 systems for which a curve 
of growth TV, b solution is available. Dotted line represents the re- 
mainder (455) for which only Lyo fits were used. The median and 
standard deviations for each sample are 27 ± 15 km (total), 
25 ± 13 km s-l (CoG), and 28 ± 16 km s"! (Lyo). Broad Lya 
Absorbers (BLAs) 6hi > 40 km s~^ represent < 20% of the ab- 
sorbers. These may indicate high gas temperatures (T > 10^ K), 
multiple components, or bulk flows within the absorber. A number 
of the broadest systems may be multiple blended components (see 
ShuU et al.(2000) and Paper 2 for discussion) so this BLA fraction 
is probab ly an upper limit. Compare with Fig. 6 in Paper 2 and 
Fig. 12 of iTripp et"an (,2008). 

absorbers is conserved, any effect on /3hi, dM /dz, or fl 
from this artifact is likely negligible and is unique to the 
H I data. 

Hydrogen is the lightest element and thus is particu- 
larly sensitive to thermal broadening. The distribution of 
doppler paramters, 6hi, provides a valuable upper limit 
on gas temperature. Figure [4] shows the distributions 
for the entire sample (650 lines) of H I measurements, 
the corrected Lya-only (455), and CoG-quahfied (195) 
subsamples. The distributions are similar for each sam- 
ple, with a median and la variation of 27 ± 15 km s~^ 
(28 ± 16 km s~^, 25 ± 13 km s~^), consistent with the 
results of Paper 2. Assuming that the entire line width is 
due to thermal broadening, 6hi — 27 km corresponds 
to a gas temperature T — 60.5 fe^j « 44, 000 K. However, 
some of the doppler width is undoubtedly due to velocity 
components and bulk flows. 

(Lehner et al.l (j2007f ) find a large number of broad Lya 
absorbers (BLAs) in many of the same sight lines with 
6hi > 40 km s~^ as measured from Lya lines. In some 
cases, they report 6hi > 100 km s~^. If the line width 
is entirely due to gas temperature, these measurements 
imply T > 10^ K. However, very broad lines are diffi- 
cult to detect against all but extremely well-characterized 
AGN continua and subject to a number of interpreta- 
tions. Only extremely strong Lya systems can be both 
broad and deep at the same time, and these are often 
not unambiguously a single component. Shallow, broad 
features are far more common, but may arise from contin- 
uum variations or instrumental artifacts such as "echelle 



ripples". We find that 15-20% of our H I systems show 
6hi > 40 km s^^, and the fraction falls off rapidly at 
higher &hi- The details of these individual detections are 
beyond the scope of this paper, but we will address these 
BLA systems in detail in an upcoming paper (Danforth 
et al., in preparation). 

3.1.2. O VI 

Our O VI results (Figure [5^,b) are compatible with 
those in Paper 1. In that work, we found a column 
density distribution with /3ovi — 2.2 ± 0.1, slightly 
steeper^ than seen in the current, larger sample (/3ovi = 
1.98 ± 0.11). The multiphase slope parameter ai4 dif- 
fers between the two samples as well (0.7 vs. 0.9). In 
Paper 1 we derived ai4 = 0.9 ± 0.1, consistent with no 
A'hi correlation; but we now find an = 0.71 ± 0.03 in 
the larger sample, consistent with some degree of corre- 
lation. However, the two studies probe different redshift 
ranges. If the current sample is limited to Zabs < 0.15, as 
in Paper 1, we have 35 absorbers with (3 — 2.26±0.20 and 
ai4 = 0.81±0.06, in closer agreement with Paper 1. This 
suggests some redshift evolution in the O VI absorbers 
between < ^abs < 0.4, as suggested in some simulations 
of the rate at which the metal-seeded WHIM increase s 
at low redshift (jPave et al.lll999l : ICen fc Ostrikeilll999l ). 

The new sample of O VI absorbers shows a strong 
power-law distribution (Figure [5Jd) . In Paper 1 , we re- 
ported hints of a turnover in the dAf /dz column-density 
distribution at logiVovi < 13.4 (VF1032 < 30 mA). How- 
ever, there were only five absorbers in these low-TV bins, 
and our conclusion was tentative. Our new study fea- 
tures 26 absorbers with logA^ovi < 13.4, and the distri- 
bution continues as a smooth power law, with no ap- 
parent turnover down to logA^ovi ~ 13-0 (W1032 > 
12.5 mA). Even the z < 0.15 subsample has six ab- 
sorbers with logA^ovi < 13.4 and also shows a similarly 
smooth power law down to logiVovi ~ 13.1, with no 
suggestion of a turnover. The continued rise in this dis- 
tribution can be used to constrain the extent of metal 
transport away from dwa rf galaxies (|Stocke et al.ll2006bt 
iTumlinson fc jPand [20051). 

As in Paper 1, we compa re the observe d O VI statistics 
to theoretical predictions. iCen fc Fand ([2006) compute 
several different sets of models of the evolution of the 
WHIM with (and without) realistic galactic superwind 
(GSW) feedback and with non-equilibrium evolution of 
major metal ions (including O VI and other WHIM trac- 
ers). They compare their models (their Figure 2) to 
the observed cumulative d/V / dz distribution at differ- 
ent O VI equivalent widths (Paper 1). They see reason- 
ably good agreement with the data in the model with 
GSW and n on-equilibrium ioniz ation. In Figure [6l we 
compare the ICen fc Fan j (|2006l ) models to the current 
dM /dz data versus equivalent width. Minimizing the 
Lya bias has generally increased dM /dz and now is best 
described by the GSW model with species in ionization 
equilibrium. The non-equilibrium model underpredicts 
O VI absorber density at all column densities. While we 
do not believe that the models are sufficiently accurate 

^ Results in Papers 1 and 2 did not use the distributed-power 
histogram technique reported here. This technique is known to 
decrease fi values by 0.05-0.10 dex, however, we have not adjusted 
the literature values. 
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Fig. 5. — Detection statistics and dST jdz for O VI (top), CIV (middle), and N V (bottom). All column densities are measured directly 
from the data (profile fits or AOD method) and not from curves of growth. 
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Fig. 6. — Observed cumulative dj\f /dz for O VI at di fferent equiv- 
alent widths, compared with several simulations by ICen fc Fane! 
(f20 06). The simulation labeled LTE (actually coUisional ionization 
equilibrium) and galactic stellar winds (GSW) matches the obser- 
vations most closely. However, note the slight excess in observed 
strong systems above W > 100 mA. This may be due to enhanced 
metallicity in the strongest systems, or it may be an artifact of the 
finite size of the simulation volume. 

to discriminate between CIE and non-equilibrium ion- 
ization, it is clear that GSW are required to produce the 
observed enhancement in number of weak, low-A'^ovi ab- 
sorbers. Changes in metallicity tend to shift the models 
up and down, and a higher metal production rate might 
match the observations more closely. There are certain 
additional systematics which must be mentioned. For 
instance, the models were calculated at z = 0, while our 
O VI observations are at (z) ^ 0.2. The WH IM fraction 
is expected to increase as a function of time (jPave et al.l 
[199 9) ; if the models were integrated over the same red- 
shift range as the observations, they would likely predict 
lower dM /dz values. Also, if we account for the full O VI 
census including lines with no measurable Lya counter- 
part, the dM /dz values would increase as discussed be- 
low. Both of these systematics would decrease the agree- 
ment between the data and models. 

We see a distribution of O VI 6-values that covers a 
wide range. This distribution is of physical importance, 
since explaining observed line width as thermal broaden- 
ing provides an upper limit to the temperature of O VI- 
bearing gas in the IGM. O VI has a thermal line width 
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(17.05 km s"^^ 



TOO 



1/2 



(11) 



where T^ax = lO''"^ K and toq = 16 to^. Turbu- 
lence and unresolved velocity components may artificially 
broaden the line profiles, but 6ovi does provide a use- 
ful upper limit to the gas temperature. For the most 
part, O VI lines were measured by a profile fitting rou- 
tine that takes instrumental resolution into account. The 
distribution looks much like those reported in Fig. 12 of 
iTripp et al.l (|2008l ) and in Fig. 6 of Paper 2. Of our 83 
O VI detections, only 11 have &ovi > 50 km s~^ and 
only four systems have 6ovi > 60 km s"^ (each of these 



four systems shows multiple components). The mean 
is (6ovi) ~ 30 ± 16 km s~^ with a median value lower 
by 3 km s~^, nearly identical to our results in Paper 2. 
Nineteen of our absorbers (23%) show 6ovi < 17 km s~^ 
(temperatures at or below the peak CIE temperature) 
while only five (6%) show 6ovi < 10 km s^"'^, correspond- 
ing to T < 10^ K, consistent with photoionization. By 
comparison. iTripp et al.l ()2008[ ) find that 37% of their in- 
tervening O VI systems have & < 10 km s~^, implying low 
temperatures in a substa ntial fraction of their absorbers. 
As shown bv lCen fc Fang. (2006.), 

one expects some frac- 
tion of the WHIM to lie at temperatures below T„iax, 
owing to time-dependent radiative cooling to T < 10^ K. 
However, Tripp et al. subdivided many of their absorp- 
tion systems into individual components, and some of 
the fits may be artificially narrow. 

As in Paper 2, we can use the distribution of H I line 
6-values to provide an upper bound on the size of IGM 
clouds. Non-virialized clouds of greater than ~ 1 Mpc 
extent would show significant Hubble expansion. Assum- 
ing that the observed line width is due to instrumental, 
thermal, and Hubble expansion effects, and that all O VI 
absorbers exist at their peak CIE temperature, we cal- 
culate 



(^) < [(feovi) 
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[i7(z)]"i - 350/ifoi kpc, (12) 



consi s tent with previous results (|Tumlinson fc Faii3 
[20051: iDanforth et all 120061) . Here, we have accounted 
for the small increase in Hubble expansion, H{z) = 
Ho [n^ (1 -f zf + f7A]i/2 w 75 km s"^ Mpc"\ at the 
mean H I absorber redshift, {z) = 0.14. 

Throughout this analysis, we assume that O VI is pre- 
dominantly thermally ionized. For photoionization to 
dominate coUisional ionization for O VI, the gas den- 
sity must be fairly low, uh ~ 10~^ cm~^. A typical 
H I absorber has A'^hi ~ lO" cm'^ and, at WHIM tem- 
peratures, a neutral fraction /hi < 10~^. This gives a 
characteristic length scale 
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1014 cm-2 
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The high ionization parameter for photoionized O VI re- 
quires extremely diffuse clouds ((5h ^ 5) and thus ex- 
tremely long pathlengths [Lh ~ 3 Mpc) for typical ob- 
served column densities. The Hubble broadening over 
this length for an unvirialized cloud, vh = ^hi ^^o ~ 
200 70 km s^i, is many times larger than the observed 
line widths of any observed absorber. Photoionized O VI 
may be feasible in high radiation environments near AGN 
at observed column densities, or in extremely low column 
density diffuse clouds, but is largely inconsistent with the 
observed IGM. 

In most instances of metal line absorption, we ob- 
serve an apparent kinematic association of H I and O VI 
(typical At; < 50 km s~i), suggesting that collision- 
ally ionized gas and photoionized gas are spatially as- 
sociated in a multipha s e system, a s simu lations predict 
(|Cen fc Ostrikeil 119991 : [Pave et al.l 1200 If ). We assume 
that the two phases have the same metallicity and dif- 
fer only in ionization state. Even at WHIM temper- 
atures, there exists a small hydrogen neutral fraction 
(/hi ~ 10~^) and thus there should exist a hot gas sig- 
nature in the Lya profile as well as our standard WHIM 
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metal ions: 
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Here, we adopt th e CIE-peak ratio /Hi/.fovi = 7.5 x 10 ^ 
at T = lO^ '^^ K (jSutherland k Dopital [1991 and scale 
A^Hi to iVovi = 10" cm-2 and Zq = 0.1 Zq. Typical 
O VI detections are A^ovi ~ 4 x IQ^^ cm~^, so the ex- 
pected hot H I line would have n'^^'^'' ~ 6 x 10^^ cm~^ 
with doppler width &hi ~ 70 km s^^, equivalent width 
Whya ~ 33 mA, and a central line depth only ~ 7% 
below the continuum (tq ~ 0.07). This is marginally 
detectable by itself in high S/N data, but it is very diffi- 
cult to separate from the stronger, narrower Lya profile 
in a multiphase absorber and might be mistaken for the 
wings of the spectrograph point spread function. Even 
a strong O VI absorber (A^'ovi — lO^"' cm~^) would 
show Wi^ya — 78 mA and a fractional depth of ^ 15% 
{to ^ 0.16). 

To investigate the existence of blended hot-plus- 
photoionized Lya profiles, we compare the 6hi distribu- 
tion of H I systems with O VI detections with those hav- 
ing clean O VI nondetections (iVovi < lO^^^'^ cm^^). For 
the 83 systems with an O VI detection, the median and 
standard deviation are &hi = 31 ± 15 km while the 
273 O VI nondetections show 6hi = 26± 13 km s~^. The 
sample means have uncertainties of ±3 km s^^ and ±2 
km s^i, respectively, which suggests that shallow, weak 
Lya lines are broadening the overall H I absorption pro- 
file. Hot H I is an ideal method of tracking the WHIM 
and does not rely on metal enrichment. However, it is 
unclear whether a difference between the two populations 
would be apparent without the O VI diagnostic. 

The O VI census in Paper 1 included only absorbers 
for which WLya > 80 mA, yet A^ovi is seen to have 
little correlation with A^hi and Iow-A'hi absorbers often 
have comparable values of Nqvi in the multiphase plot. 
What effect does this "Lya bias" have on the WHIM 
census? Weak Lya absorbers {Whya < 80 mA) make 
up 45% of our current survey, while 19 out of 83 O VI 
detections (23%) were seen in weak Lya lines, which we 
would not have considered in Paper 1. In this survey, 
for the FUSE sample of 35 absorbers at z < 0.15, we 
find dJ\f /dz = l^tl down to 30 mA and dN /dz = eot^g 
down to 10 mA. These compare to our finding in Paper 1, 
dN /dz = 17 ± 3 (down to 30 mA) where weak H I sys- 
tems were not taken into account. Our estimate, ilovi 
(Table 5) is also ^ 25% higher in the low-z sample than 
was seen in Paper 1. 

Our current survey does not fully correct for Lya bias, 
since we still require an H I detection at some level 
to trigger a search for O VI. I Tripp et"al] ()2008[ ) per- 
formed an ambitious "blind" O VI survey of 16 AGN 
sight lines including many of those in this work. They 
find dAf /dz = 18.3^2 g down to 30 mA, compared to our 
dM /dz — 15.0^2 Q to the same limit (z < 0.4 sample). In 
the same sub-sample, we also find dj\f /dz = 40^g^ down 
to 10 mA {dJ\f/dX— 54jlg^). Ignoring systematic differ- 
ences between the two surveys, such as our respective 



criteria for absorber path length (Az) and S/N-limited 
equivalent width limits, this implies that an additional 
~ 20% of O VI absorbers do not have any measurable 
corresponding H I absorption. This effect was already 
apparent in a qualitative sense in the Iow-A'hi end of the 
multiphase plot (Figure lb of Paper 1), where we found 
many absorbers with Nqvi > Nm. When combined with 
the ^ 20% effect implied by the differences between this 
work and Paper 1 noted above, we conclude that "Lya- 
bias" can have a significant effect on a complete WHIM 
baryon survey. 

Following the example of iTumlinson fc Faiij (|2005f l , 
we can use the O VI line frequencies, dM /dz, to con- 
strain the distribution of metals in the low-z IGM. Re- 
ferring to their Figure 2, we find that our value for the 
z < 0.4 sample, dM/dz — 40lg*, integrated down to 
10 mA, is consistent with a mean spatial extent of met- 
als out to 250-300 kpc from galaxies in the Sloan Dig- 
ital Sky Survey (SPS S). In this picture, and that pro- 
posed bv lStocke et al.l ()2006b), metal enrichment occurs 
around dwarf galaxies of SDSS r-band absolute magni- 
tudes down to Mr- « -17 to -18 (0.03-0.05 L*). 

3.1.3. CIVandNV 

The ions N V and C IV are both Li-like, with 
strong UV resonance (2s — 2p) doublets. Like O VI, 
they may provide good WHIM tracers, at longer wave- 
lengths accessible in the HST spectral band. These ions 
have be en studied toward many G alactic high-velocity 
cloud s (ISembach fc Savagd 1 1994 llndebetouw fc Shulll 
I200I ISembach et al.l 120031 : IColhns et all 120031 12007^ 
However, to date, no large, low-redshift IGM surveys 
have been conducted in these ions. C IV is a strong 
absorber of an abundant element, but the ST1S/E140M 
data arc limited to Zabs ^ 0.116, at which redshift C IV 
A1548 shifts out of the passband. The longer wavelength 
STIS/E230M grating covers 0.01 < zciv < 1-0, but this 
band was not often used in AGN observations. 

In our new survey, we detected 24 C IV absorbers 
in the ST1S/E140M data over a total redshift path- 
length Az = 2.42 and inferred {dN /dz)ciw — IQtl and 
rJciv = (7.67 ± 1.45) x 10-8 down to > 30 mA. At 
high redshift (1.5 < z < 3.1) IScannapieco et all ()2006[ ) 
inferred ilciv = 7.54 x lO"* wit h no sign of evolution over 
that redshift range. Similarly, ISongailal (jl997, 20 0l|) in- 
ferred f^civ ~ (2-7) x 10"^ at z = 3.0-3.5. iPettiniet all 
12003) reported l^civ = (4.3 ±2.5) x 10"* at z ~ 5 using 
similar column density limits as are used in this paper. 
The consistency of our low-z measurement to those in 
the early universe implies little or no evolution in C IV 
since reionization. This may imply a changing metal en- 
richment rate and/or changing ionizing fiux with red- 
shift. However, since the ionization corrections depend 
on both photon flux and gas density, further speculation 
is beyond the scope of this paper. For carbon metal- 
licity Zc = 0.1 Zq, where (C/H)q = 2.45 x 10-^ and 
/civ ~ 0.29, the maximum ion fraction at logT^ax = 5.0 
in GIF, r^civ corresponds to a total baryon fraction asso- 
ciated with G IV of r^lcM^ = (2.6±0.5) X lO'^ or 25-30% 
that seen with O VI. 

As with O VI and H I, the G IV absorbers fol- 
low a power-law distribution in column density with 
/5civ = 1.79 ± 0.17 (Figure [5j:,d), consistent with litera- 
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ture values at high redshift, civ = 1-5 — 1.8 (jSongailal 
119971 120011 : lEmson et al1l2000l ). This slope is shallower 
than is seen in O VI and more consistent with those of 
lower-ionization species. However, the uncertainty on 
Pew precludes any definite conclusions as to whether 
C IV exhibits a steep, WHIM-like distribution in col- 
umn density, or a more moderate, H I-like distribution. 
The ionization energy required to produce C IV is only 
3.52 ryd (vs. 8.37 ryd to produce O VI), thus C IV 
may be produced partl y by photoionization by AGN 
and hot stars at low z (jGiroux fc Shul]||l997T ). On the 
other hand, we see little correlation between A^civ and 
the column densities of lower ions, even C III (see be- 
low). In contrast, there is some degree of Nqiv-Novi 
correlation, suggesting that C IV may primarily trace 
shock-heated material. The baryon fraction traced by 
C IV is less than half that of O VI for a given equiv- 
alent width limit, consistent with the solar abundance 
ratio, (G/0)p, = 0.5 ± 0.07 (|Allende Prieto et al.l[200l 
lAsplund et al.l [2005( 1 . and the implied sub-solar (C/0) 
ratio seen in the IGM (Paper 2). In the pure thermal 
(collisional ionization) interpretation of C IV, the cooler 
portions of WHIM (at T ~ 10^ K) may contain a smaller 
portion of the baryons than at T '--^ 10^'^ K, consistent 
with more efficient cooling at temperatures near or below 
the peak of the cooling curve. 

The rest wavelengths of the N V doublet (1238.821, 
1242.804 A) are well placed in the STIS/E140M data 
for redshift coverage out to Zabs ~ 0.4. However, the 
solar abundance of nitrogen is less than 20% that of oxy- 
gen. Moreover, nitrogen is a product of secondary nu- 
cleosynthesis, and its abundance relative t o C and O is 
observed to be sub-solar in Gal actic HVCs ([Gibson et al.l 
[200HlCollins et al.l[2003l [20f)l . We detected 24 N V ab- 
sorbers in one or both lines of the doublet (Figure [5^,f ) 
over a total pathlength Az — 5.30. Because of the lower 
relative abundance of nitrogen compared to carbon or 
oxygen, the detections are never strong, and we derive 
dAf/dz = 2 ± 1 down to 30 mA. Weaker N V lines are 
more numerous, and our total sample has d/V /dz = 7^.1 
down to 10 mA. We detected N V and O VI together 
in 11 systems, and their column densities are well cor- 
related. The ionization potential to produce N V is 
5.69 ryd, and its peak CIE ion abundance occurs at 
logTmax = 5.25, near to tha t (logT^ax = 5.45) for O VI 
([Sutherland fc Dopital Il993f ). We also note that these 
high ions are not necessarily in equilibrium, and non- 
equilibrium ionization effects may alter /ion- We infer a 
value of il|Qj2/^ft intermediate between those of O VI 
and C IV at both 30 and 10 mA (Table [T2|). We con- 
clude that N V is a reliable tracer of WHIM gas based 
on its correlation with O VI and steep power-law slope, 
although with sizable error bars. Future studies of these 
ions with the Cosmic Origins Spectrograph on Hubble 
have great promise, not only to trace the WHIM but to 
derive nucleosynthetic signatures from a range of heavy 
elements (e.g., C, N, O, Si, Fe). 

3.1.4. cm 

We made 39 C III detections, yielding dN/dz^ lOtl, 
{dM/dX = 12^2) down to 30 mA. The column den- 
sity distribution shows a power-law slope of /3ciii = 
1.79 ±0.10 (Figures [7^, b). Absorber frequency is com- 



patible with that found in Paper 2, but spectral slope 
appears higher by ~ 0.2 dex. A subsample of absorbers 
at Zabs < 0.21, comparable to the sample used in Pa- 
per 2, shows 28 absorbers with a somewhat higher ab- 
sorber frequency, dN jdz^ 14^3 (dM jdX = 14^4), and a 
steeper column density distribution, /3ciii = 2.06 ± 0.15. 
The reason for this discrepancy may be related to our 
inclusion of weaker Lya systems, which were not studied 
in our previous work. Since much of the C III is ex- 
pected to arise in the same photoionized gas traced by 
H I, we would expect that the inclusion of weaker H I 
systems in the survey would produce a larger number of 
C III systems as well. However, this does not explain the 
steepening of the column-density distribution seen in the 
low-z subsample. 

In Paper 2, we saw a reasonably strong correlation be- 
tween A'hi and A^'cm. The multiphase ratio was fitted 
with a line of slope a 14 = 0.73 ± 0.08 and scale factor 
Ci4 — 1.06 ± 0.04. From this, we incorrectly derived 
a metallicity for intergalactic carbon of Zc = 0.12 Z0, 
in good agreement with Zq — 0.09 Zq derived in Pa- 
per 1. The correct derivation was lower by a factor of 
ten and should read Zc = 0.012 Z©, for a peak CIE 
ionization fraction /cm = 0.83. This inferred metal- 
licity would increase if we chose a lower value of /cm 
more typical of photoionized gas. While the corrected 
value is substantially different from the derived oxygen 
metallicity in the IGM, the implied abundance ratio, 
(C/0)iGM ~ 0.1 (C/O)©, is compatible with the ob- 
served C III/O VI ratio (see § 4.4 and Figure 11 of Pa- 
per 2). With more C IV data, we could use both ions, 
C III and C IV, as a joint constraint on ionization cor- 
rections and carbon abundance. 

The behavior of C III in our new study is consistent 
with the (corrected) previous results in Paper 2. The 
correlation of iVni and A^cm is still strong, and metal- 
licity implied by the multiphase fit {Zq/Zq w 1.5%) 
is similar to previous results for both the full sam- 
ple of 39 absorbers and the 28-absorber low-z subsam- 
ple. When adjusted for the mean absorber redshift of 
each sample and corrected by a range of ionization frac- 
tions, /cm =0.2-0.8, for either photoionization or col- 
lisional ionization equilibrium, these values imply Zq = 
(0.015 ± 0.002) ZQ(/cm/0.83)-i or Zq/Zq = 1.5-6%. 

3.1.5. Si III, Si IV, and Fe III 

Intergalactic Si III, Si IV, and Fe III at low redshift 
have received little or no attention. Along with C III dis- 
cussed above, all three ions are are presumed to trace pri- 
marily metal-enriched, photoionized material. All have 
peak CIE temperatures below 10^ K and ionization po- 
tentials > 3 ryd. The cooling times at these tempera- 
tures are relatively short, and even at the low densities 
typical in the IGM, any coUisionally heated material at 
T < 10^ K should quickly cool. Such transient situations 
are likely places for non-equilibrium ionization effects. 

Because Si III A 1206. 500 is a very strong transition, 
STIS observations are sensitive to extremely low column 
densities. Despite the lower abundance of silicon rela- 
tive to more common metals ([Si/O]© ~ —1.4), Si III 
is expected to be the dominant ionization stage in the 
photoionized or coUisionally ionized IGM over a broad 
range of temperatures and photon energies. We there- 
fore expect it to be an effective tracer of IGM enrichment. 
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Fig. 7. — Same as Figure[3]for C III (top), Silll (second row), SilV (third row), and Felll (bottom). 
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complementary to C III or O VI. We detect Si III in 53 
systems (Figures [7b,d) for dAf/dz= Qt\ [dAf /dX = 7t\) 
down to 30 mA or dN /dz= Ut^ [dN / dX = IQtt) down 
to 10 mA, about 40% of the value for O VI at either limit. 
The column-density distribution, {d'^Af/dz dNsun), fol- 
lows a clear power law with /3siiv = 1-80 ± 0.09, consis- 
tent with Pciu and /3ciVj slightly higher than /3hi- 
We made 20 Si IV detections and derived dJ\f /dz= 4^^ 

(dAf/dX = 5ti) down to 30 mA (Figures [7l3,f). The 
column density A^siiv does not show a particular corre- 
lation with A'hi, as do A^ciii and A^siiii- The column 
density distribution follows a rough power law and a fit 
to the cumulative distribution gives /3siiv = 1-92 ± 0.17, 
higher than other low-ionization species. We measure 
f]giiv = (8.3 ± 1.3) X 10"^ down to 30 mA. Fo r absorbers 
at at high redshift, 'Scannapiec o et all (j2006l ) measured 
fisiiv = 6.0 X 10"^ to logiVsiiv> 12 at (z) = 2.4, a fac- 
tor of 14 lower than our most sensitive limit at (z) ^0.1. 
They see no evidence for evolution in Si IV in the redshift 
range 1.5 < 2 < 3.1. Similarly, ISongailal (|1997l l2001h 
measured flsnv = (7 - 30) x 10~^ at 1.5 < z < 5.5. 
This factor of 2.5 — 12 increase in the Si IV fraction be- 
tween the early and modern universe is opposite to the 
behavior of C IV above, which remains essentially un- 
changed. This may imply some combination of changes 
in the metagalactic ionization field and/or IGM enrich- 
ment. 

Fe III is a relatively weak transition (/ = 0.0544, 
iMortonI |2003[ ) of a relatively low- abundance element 
([Fe/O]0 « —1.4) and thus it comes as no surprise 
that we make only 14 detections (Figure [7^,h) with 
dN/dz down to 30 mA. The power-law slope is the 
steepest of any in our survey (/3fcIII = 2.2 ± 0.4) but 
with such large uncertainty as to render any compari- 
son speculative. As with iVovi, -^FoIII shows only cur- 
sory correlation with A^hi and covers a range of 1.5 dex 
in column density, while A'hi is seen over almost three 
decades. This evidence is exactly what we used in Pa- 
per 1 to argue that O VI and H I were tracing different 
IGM phases. However, given the low ionization energy 
required to produce Fe III (16.18 eV), it is likely that 
Fe III and H I arise in similar photoionized gas. We may 
be seeing preferential iron enrichment in material that 
has been shocked and cooled to ambient temperatures. 
We will return to these issues and their redshift depen- 
dence in a future paper that deals with IGM absorber 
evolution. 

3.2. Ion Ratios and Comparison with Models 

With such a rich data set, many absorbers show de- 
tections in more than one ion, allowing us to compare 
observed line ratios with predictions from a set of sim- 
ple CLOUDY models. In Paper 2 we compared the ob- 
served C III/O VI ratios for 13 absorbers with models 
and concluded that a single thermal-plus-photoionized 
phase could not realistically account for the observations. 
We found that C III and O VI absorbers probably in- 
cluded both photoionized and coUisionally ionized gas, 
although O VI was likely to be mostly coUisionally ion- 
ized. 

In this work, we use the same set of models: a 
four-dimensional grid calculated using CLOUDY v96.01 
with parameters temperature T, ionization parameter 



U , metallicity Z (relative to solar values), and hydro- 
gen number density n (cm~'^). We simulate the effects 
of coUisional ionization by setting the cloud to a con- 
stant temperature in each model. Low temperatures, 
T K, (1 — 3) X 10** K, approximate a pure photoion- 
ization model, while low photoionization parameters, 
log U < —4, are equivalent to a pure CIE model. These 
models were calculated for a 400 kpc thick optically thin 
slab of gas illuminated on one side by an AGN radia- 
tion field with a power-law continuum F^, cx in the 
EUV ( 1-4 ryd) and soft X-ray (10-22 ryd) with a, = 1.8 
(Haar dt fc Madaul 119961 : iFardal et all 119981 : iTelfer et all 
12002) . Since a more realistic isotropic radiation field il- 
luminates the slab from all sides, the actual photon den- 
sity is a factor of 4 higher (-1-0.6 dex) than that output 
by the models; we have adjusted all citations to refiect 
this. We also computed a second grid of CLOUDY mod- 
els assuming n = 10~^ cm~^, Z = 0.1 Zq, and pure 
photoionization (no thermal ionization component), but 
we allowed the ionization parameter U and spectral in- 
dex as to vary b etween 1 < < 3 (jTelfer et al.ll2002l : 
iShull et al.l[200l . The equilibr ium temperatures of the 
pure photoionization models was logT = 4.1 — 4.3. 

Observed Si III and Si IV show a strong correlation in 
column density {R = 0.85, A^siiii oc A^siiv^ °=^° \ Fig.[li), 
suggesting that the two ions trace the same gas. Since 
Si Ill/Si IV are adjacent ionization states of the same 
element, we can avoid complications arising from abun- 
dance and metallicity and use A^siiii/A^Siiv as a tracer 
of the ionization state of the cloud. The observed ra- 
tios show some scatter, (log (A^sim/A^siiv)) = 0.1 ± 0.4, 
but all 12 systems in which both ions appear can be 
explained as a single-phase absorber, either collision- 
ally ionized at T = lO'^-^^O-^ K or photoionized with 
\ogU = — 1.6± 0.7. Variations in model metallicity and 
density do not produce significant shifts in allowed pa- 
rameter space, nor does changing the photoionizing spec- 
tral index as significantly alter the (U, T) solution. Un- 
fortunately, the doppler &- values for the lines are of lit- 
tle help in constraining the temperature solution; the 
median line width for these Si ions is ^ 30 km s~^, 
corresponding to T > 10^ K for a case of pure ther- 
mal broadening. We conclude from this that the low- 
ionization material is consistent with photoionized clouds 
with \ogU ~ —1.6. Similar CLOUDY models and silicon 
ion ratios in Galactic high velocity clouds (HVCs) show 
log { / ^ — 3, presumably due to much higher gas dens ities 
(see lCollins et~all[200a [200^ iWestbrook etall[200l . 

C III and C IV are the other pair of adjacent ions 
in this work. Detections in both species occur in only 
eight absorbers (Figure [8)3) , and the column densities 
are not well correlated [R ^ 0.5). Still, the detections 
are consistent with a 1:1 column-density ratio, and the 
absolute scatter in log (Acm/A'civ) is actually smaller 
than for log (A'siiii/Asiiv)- Compared to the models, the 
range of observed N cm/ New is consistent with either 
logC/ = -1.4 ± 0.5 or logT = 5.0 ± 0.1, near the peak 
CIE fractional abundance of C IV. However, since A'civ 
and A^ciii do not show the same degree of correlation 
present in the Si ions and since we believe C IV to be 
coUisionally ionized in at least some cases, we hesitate 
to assign too much significance to the C III/C IV model 
comparisons. 

Because of the expected redshift evolution of the 
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Fig. 8. — Detections of both Si III and Si IV (left) and CIII and CIV (right). The column densities of the Si ions are well-correlated 
{R = 0.87), and the linear fit (dashed) has a slope l.litO.l, consistent with a 1:1 ratio (dotted). The correlation is much less for the carbon 
ions (R = 0.52), although the data are also consistent with a 1:1 ratio (dotted). The absolute scatter for each pair of ions is similar. 
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Fig. 9. — Line ratio plots for the same ionization stages of different elements. Six absorbers with both CIV and Si IV detections show 
some correlation (left), but the general trend is consistent with Nciv oc A^silV (dotted). Si III and Felll (right) appear in seven absorbers 
and are extremely well-correlated over several decades, but with iVpejjj oc ^gjfn (dashed), not 1:1 proportionality (dotted). 



thermal state, ionization state, and metallicity of the 
IGM, it is useful to compare our low-z results with 
those seen in the early universe. High redshift {z > 
2) Si IV is often compared to C IV as a probe of 
the m etagalac tic ionizing b a ckgrou n d ([S ongaila' 119981: 
Kim Cristiani. fc D'Odoricol 120021 : IScannapieco et al.1 
20061) , or as a measur e of the a-proce ss versus RGB 
stellar nucleosynthesis (|Mc WilliamI 1 1 997i ) . however both 
analyses assume that both ions are photoionized. For 
young nucleosynthetic systems, one expects to observe 
abundance enhancements in a-process elements (Si, O) 
relative to (C, Fe). However, deriving accurate (Si/0) 
and (Si/C) abundances is difficult, owing to the de- 
pendence of ioni zation corrections on the shape of the 
ionizing spectra (jGiroux & ShuU [19971 ). We see six ab- 



sorbers with both Si IV and C IV detections with — 1.3 < 
log (A^Siiv/A^civ) < —0.3, similar to the range seen at 
higher redshift. The correlation is moderate {R = 0.6), 
but the general trend is consistent with a linear relation- 
ship between the two species (Figure UK)- 

The ion pairs Fe Ill/Si III and C Ill/Si IV have similar 
ionization potentials and thus their column density ratios 
should be relatively insensitive to the ionizing spectrum. 
These two pairs provide probes of relative C/O/Si/Fe 
abundance ratios in the IGM. C III and Si IV are seen 
together in six absorbers and are poorly correlated. How- 
ever, the observed (log(A^ciii/-^Siiv)) = 0.7 ± 0.6 is 
slightly lower than model predictions at \ogU sa —1.5, 
suggesting, if the single-phase assumption is used, an 
overabundance of silicon/carbon by ^ 0.3 dex with re- 



22 



spect to solar values. 

Fe III and Si III are seen together in 7 systems 
and their column densities are well-correlated {R > 
0.9) though with iYpeiii oc Afgifif"-^ (Fig. [Hb). The 
modeled iVFoiii/-/Vgini values depend on U, but for 
\ogU w —1.5 and solar [Fe/Si], the observed values, 
(log(iVFciii/-^Siiii)) = 0.5 ± 0.5, are too high by a factor 
of 2-3 dex. If Fe III and Si III truly arise in the same gas, 
normal photoionization corrections require an enormous 
overabundance of iron in the IGM, [Fc/Si]iGM = 2.6±1.2. 
A similar iron overabundance is found by comparing 
Fe III/C III. 

The basic issue here arises because of the low ioniza- 
tion threshold (30.63 eV) and high cross section for pho- 
toionizing Fe III to Fe IV. If Fe/Si and Fe/C are in so- 
lar abundance ratios, it is difficult to obtain iVpoiii > 
^Siiii- In many of the IGM absorbers detected in both 
ions (Fig. [SJj), the observed column density ratio is 
[Fe Ill/Si III]= 0.5 ± 0.3, whereas CLOUDY models pre- 
dict much lower ratios. The solution to the Fe III prob- 
lem may lie in details of the assumed AGN spectrum; 
we use a continuous power-law spectrum from 1-22 ryd, 
which leads to Fe III being photoionized to Fe IV. A 
power law with a slope break at 3-5 ryd, as proposed by 
iHaardt fc Madaul ()2001[ ) , results in relatively less flux at 
higher energies and boosts the predicted Fe III columns. 
However, this would make the production of photoion- 
ized C IV and other high ions even more difficult. With 
more data and a larger survey, we might be able to dis- 
tinguish among several competing effects: metallicity en- 
hancements of [Fe/Si] due to nucleosynthesis, changes in 
the ionizing radiation field changes around 3-5 ryd, and 
the relative portions of Fe III, Si III, and C III in various 
ionization phases. 

The CLOUDY models can also be applied to the high 
ions, believed to reside in the WHIM. The correlation is 
much better for the eleven systems with both N V and 
O VI absorption (R = 0.8, Fig.[lOK). The observed best 
fit has a slope of 0.75 ± 0.08, but a fit with slope of 1 is 
quite plausible with an offset of — 0.74±0.03 dex, compa- 
rable to the solar abundance, log(A^/0)Q ~ —0.7. Un- 
fortunately, these three ions are neither of the same ele- 
ment nor do they have overlapping ionization potentials, 
so there is an abundance/ionization degeneracy in the 
{U,T) solutions. Assuming a solar N/0 abundance, we 
find that the observations are consistent with collision- 
ally ionized systems at T = lo^-ssio io K, midway be- 
tween the peak GIF abundance temperatures of the two 
ions. Changing the N/0 abun dance to the sub-solar val- 
ues typically seen in HVC gas (iGibson et al.l[2000L [200ll : 
[Collins ct al. 2003, 200"^ does not change the tempera- 
ture solution appreciably. 

Observed N V/0 VI ratios are also consistent with 
pure photoionization models with \ogU = — 0.4 ± 0.5. 
However, this is a factor of ^ 10 higher than that de- 
rived from the low ion ratios. This pattern is familiar 
from the O VI/C III results in Paper 2, where G III sug- 
gested logf ~ —1.5, while photoionized O VI required 
\ogU ^ —0.5. Altering the relative N/0 abundance 
in the IGM varies the ionization parameter solution. 
Matching the observed ratios to the derived logU = — 1.5 
requires a relative N/0 overabundance of 0.3-1.3 dex, op- 
posite to the trend observed in HVGs. Varying the slope 



of the ionizing spectrum will also alter the predicted ra- 
tio, but an extremely hard incident spectrum (as < 1) 
is required for a \ogU = —1.5 solution. Based on ob- 
servations of nearby AGN, a spectral slope this steep is 
deemed unrealistic. 

C IV and O VI appear together in eleven systems 
though the correlation is not as clear as with N V vs. 
O VI (i? = 0.6, Fig. Wp)- The best fit slope is 0.7 ± 0.1 
though a 1:1 fit with an offset of —0.64 ± 0.04 is plausi- 
ble. This offset may be a manifestation of the sub-solar 
C/0 ratio implied in Paper 2, or it may reflect the faster 
cooling times experienced by lower temperature shocked 
gases. 

The model ([/, T) solution is similar to that for 

N V/0 VI above; log[/ 0.5 or T - io5-4±0-i k for 

a solar C/0 ratio. Again, matching observed C/0 to 
a logC/ ~ — 1.5 solution requires G to be overabundant 
with respect to O by 0.6 dex, opposite the subsolar 
G/0 ratio implied in Paper 2. 

3.3. WHIM Cosmology and Metallicity of the IGM 

Integrated down to the lowest level of detectable ab- 
sorbers (W ^ 10 mA), our survey of three high metal 
ions estimates that they trace baryon contents of approx- 
imately « 9% (O VI), 5% (N V), and 3% (C IV) 
of the total. These estimates scale inversely with an as- 
sumed IGM metallicity of 10% solar and the adopted ion- 
ization corrections characteristic of maximum abundance 
in GIE: /ion = 0.22 (O VI), 0.24 (N V), and 0.29 (C IV). 

We see that the WHIM baryon fraction, f^jcM' implied 
by O VI absorption is greater than that from N V, which 
in turn is greater than that from G IV. However, these 
differences should not be taken too seriously, since there 
are many assumptions that go into these estimates. First, 
the ionization fractions may differ from GIF, particularly 
if photoionization plays some role in the C IV absorbers. 
Second, these ions have different sensitivity limits in col- 
umn density owing to line strengths and their location 
at different wavelengths. Alternatively, the differences 
may represent a true gradient in fiwHiM as a function of 
temperature. The latter interpretation is plausible since 
O VI, N V, and G IV have peak GIF temperatures of 
3, 2, and 1 x 10^ K, respectively, and radiative cooling 
rates h ave a steep dependence on temperature in that 
regime ([Sutherland &: Dopitalll993l ). Given its poor cor- 
relation with O VI and N V, many C IV absorbers may 
be photoionized gas rather than thermal in nature. Thus, 

^WHiM should probably be taken as an upper limit un- 
til we can refine the multiphase nature of the absorbers 
with additional data from larger C III and C IV surveys. 

At our most sensitive O VI limit {W > 10 mA), we 
find ilvvHiM = 8.6 ± 0.8%. This is higher than the value 
quoted in Paper 1 (^ 5%), probably because we now in- 
clude weaker Lya systems (the Lya-bias) and have better 
statistics in the weak O VI absorbers themselves. How- 
ever, the power-law slope of the WHIM ions is steep, 
with a similar amount of mass in the numerous weak 
absorbers as in the infrequent strong absorbers. Future 
WHIM surveys sensitive to lower A'^ovi limits may well 
increase the baryon mass fraction in th e WHIM closer to 
the 30-50% pred icted by simulations (jDave et al.lll999t 
iCen fc Fang|l2006ll . 

X-ray probes of the WHIM at T > 10^ K, where 
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Fig. 10. — Correlations between highly ionized species. The column densities of NV and O VI (left) are well correlated in eleven absorbers 
(i? = 0.81), and the best fit (dashed) has a slope of 0.75 it 0.08. However, a 1:1 fit is plausible (dotted), with an offset of —0.74 it 0.03 
dex, consistent with log(N/O)0 ~ —0.7. CIV and OVI appear together in eleven absorbers (right), but the correlation is only moderate 
{K = 0.62). Much like NV vs. OVI, the best fit slope is ~ 0.7 (dashed) but the data are are consistent with a 1:1 relationship with an 
offset of ~ -0.6 (dotted). 



a considerable fraction of WHIM is predicted to exist, 
are not yet available with the same fidelity as FUV 
WHIM tracers (O VI, etc). For example, O VII A21.60 
is 9.1 times weaker (/A) and O VIII A18.97 is 17.4 
times weaker than O VI A1032. Partially offsetting 
these atomic ratios are the CIE peak ionization frac- 
tions of O VII and O VIII, which are 3-4 times that 
of O VI. However, the primary contribution to the 
low X-ray spectral sensitivity comes from the low res- 
olution of the grating spectrographs aboard Chandra 
(700 km s-i at 20 A) and XMM-Newton (-1000 km s~'^), 
compared to those on Hubble (7-20 km s~^) or FUSE 
(20 km s~^). Thus, the O VI surveys probe hot gas 
down to A^ovi ~ 10^'^ cm~^, whereas the recent claimed 
O VII and O VIII detections are at the leve l iVoviii ~ 
5 X 10^^ cm-^ (toward P KS 2155-304, iFang et all 
I2002t IFang. Canizares. fc Yao 2007 ) and iVovii « 1 x 
10^^ cm-2 (toward Mrk421. iNicastro et all 12005') with 
significant uncertainties. However, these redshifted IGM 
absorbe rs have not been seen in XMM/Newton o bser- 
vations asmussen et al.ll2007l:lKaastra_e t al."2006') and 
their verification is therefore in doubt ([Bregman 2007). 
Thus, for the moment, the X-ray detection statistics of 
remain highly uncertain while FUV spectroscopy pro- 
vides our best (only?) measurements of the WHIM, lim- 
ited to the regime with T < 10^ K. 

In the photoionized H I gas in the Lya forest, 

iPenton. Stocke. fc ShuU ((200l found that 29 ± 4% of 
the local baryon mass could be accounted for down to 
log TVhi = 12.5 (Wlyq = 17 mA). We determine values of 
ri for H I via equations[51and[TUland present the results in 
Table [T3l and Figure [TT] for different logiVni values. Our 

survey of 650 Lya absorbers gives n[f^/nb = 29 ± 4% 
integrated down to the same column density limit as Pen- 
ton et al. but with ^ 3 times higher pathlength Az. 
Both results are consistent with cosmological simula- 
tions, which find that n'^^/ilb = 20-40% ([Pa 
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s et al.l 



fT999l[200ll: iCen fc Ostrikal[T99l[Ce"n fc Fang|[200l . 

Low ions such as C III and Si III likely trace the sub- 
set of the Lya forest that is enriched by star forma- 
tion and galaxy feedback. Using these ions as probes 
of metallicity will allow us to look for evolution of metal- 
licity with z and iVni. We can also look for patterns 
indicative of nucleosynthesis near galaxies and effects 
on ion state arising from multiphase photoionized and 
collisionally ionized gas. In the doubly-ionized species, 

we see n'^^'^/nb = 1-4% (Table [HI), well below the 
— 30% value for the Lya forest as a whole if we take 
the fiducial CIE values for /ion listed in Table 3. The 
assumed ion fractions for C III, Si III, and Fe III are 
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Fig. 11. — Differential Hhi as a function of logA'^jjl for 650 Lya 
absorbers calculated by the methods of Penton et al. (2000, solid) 
and Schaye (2001, dotted). The logAf < 12.5 bin contains only 
17 absorbers and is not included in our calculations of cumulative 
Urn in Table [T3l 
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TABLE 13 

Baryon Content of the Local Lyo Forest 



log A^Hi Range 


M 








f2Lya/f^6 (%) 


12.5 - 13.5 
13.5 - 14.5 
14.5 - 15.5 
15.5 - 16.5 


373 
206 
50 

4 


0.0053 ± 0.0015 
0.0042 ± 0.0004 
0.0031 ± 0.0006 
0.0004 ± 0.0002 


11.7 ±3.3 
9.3 ±0.8 
6.9 ± 1.3 
0.8 ± 0.4 


0.0062 ± 0.0019 
0.0033 ± 0.0003 
0.0017 ± 0.0003 
0.0002 ± 0.0001 


13.7 ±4.1 
7.3 ±0.7 
3.7±0.7 
0.3 ±0.2 



12.5 - 


14.5 


579 


0.0096 ± 0.0016 


21.0 ±3.4 


0.0096 ± 0.0019 


21.0 ±4.2 


14.5 - 


16.5 


54 


0.0035 ± 0.0006 


7.7 ± 1.4 


0.0018 ± 0.0003 


4.0 ±0.7 


12.5 - 


16.5 


633 


0.0131 ± 0.0017 


28.7 ±3.7 


0.0114 ±0.0019 


25.0 ±4.2 



Method of Penton, ShuU, & Stocke (2000)*= Method of Schaye (2001) 



all close to unity. Lower values of /ion, as might be 
the case in non-equilibrium ionization, would produce 

smaller baryon fractions and we conclude that J^Iq^ is 
an upper limit in the case of the doubly- ionized species. 
Likewise, the baryon fraction probed by Si IV at low z 

is surprisingly high, illgJ^-'/rib = 8± 1%, comparable to 
that of O VI. However, Si IV is most likely photoionized, 
and thus the CIE ion fraction /siiv ~ 0.2 assumed in 
this calculation may be too low. These two limits sug- 
gest 1% < {n^^^/nb) < 8%, or that roughly 3-25% of 
the local Lya forest by mass is enriched at a currently 
detectable level. More sophisticated models may shed 
light on this problem. 

As in Papers 1 and 2, we estimate the metallicity of 
IGM material based on multiphase ratios (Table [TT|) . 
The values generated cover a wide range from Z-pc — 
0.7 Zq to Zcm = 0.015 Zq. With Zqvi and Zcm we 
reproduce the result, [C/0] « 0.1 [C/0]q, reported in 
Paper 2, however, C IV reports Zc = 0.055, a factor of 
3 higher than estimated by C III. Similarly, adjacent Si 
ions disagree on IGM Si abundance by a factor of eight. 
As with ^^Iqm above, each metallicity estimate scales as 
[/ion]^^, so that the inferred metallicities for the low ions 
may be upper limits. However, this puts the adjacent-ion 
metallicity estimates even further from agreement. With 
larger data bases, it may be better to perform a joint 
analysis of ion pairs such as C III/C IV and Si Ill/Si IV, 
to arrive at consistent values for log U and Z. 

If we restrict ourselves to species with the most detec- 
tions (O VI, C III, Si III) we see that the WHIM tracer 
(O VI) shows Z « 0.15 Z© while the Lya forest tracers 
(C III, Si III) have a lower metalhcity {Z ~ 0.02 Zq). 
The pattern can be forced a bit by including C IV 
{Z ~ 0.05 Zq). This would imply that shocked material 
traced by high ions tends to have somewhat higher metal- 
licity than the ambient IGM, as one would expect if both 
shocks and enrichment were products of galactic winds 
and star formation feedback. However, N V, Si IV, and 
Fe III break the metallicity-ionization potential pattern 
with anomalously hi gh values. Future nearest-neighbor- 
galaxy studies fe.g.. |Prochaska et al.ll2006t IStocke et al.l 
I2006bh may help make sense of these confusing mul- 
tiphase ratio metallicity estimates. In particular, nu- 
cleosynthetic effects may show up in galactic outflows 
which may be partially responsible for IGM shock heat- 
ing traced by O VI. 

4. CONCLUSIONS AND SUMMARY 



We present the results of the largest low-redshift IGM 
survey to date. In total, we analyzed 650 H I absorbers at 
z < 0.4 along 28 AGN sight lines. For each IGM system, 
we measured detections or upper limits in 13 transitions 
of 7 metal ions as well as H I (Lya, Ly/3, O VI, N V, 
C IV, C III, Si IV, Si III, and Fe III). 

Our O VI results reinforce Danforth & ShuU (2005) re- 
sults but over a broader redshift range, and to a slightly 
deeper sensitivity limit. We found 40 O VI absorbers in 
Paper 1 at z < 0.15 and expand the sample to 83 O VI 
absorbers at z < 0.4 in this work. The dAf /dz turnover 
below log A^ovi = 13-4 seen in Paper 1 is not reproduced 
in the larger dataset, at least down to logTVovi = 13-0. 
Detection statistics from the low-z bin of this survey are 
consistent with equivalent values from Paper 1. Compar- 
ison of dhf / dz, ai4, and 17 between high and low redshift 
bins hits at z-evolution in dM jdz, ai4, and derivative 
quantities. 

The N V detections are not as common as O VI (24 
vs. 83) owing to the lower cosmic abundance of nitrogen. 
However, N V shows excellent correlation with O VI in 
column density and displays many of the same charac- 
teristics (steep power-law slope, lack of correlation with 
H I) and we conclude that N V is a reliable tracer of 
WHIM material. C IV is a strong transition of an abun- 
dant metal, but our detection statistics are hampered 
by the limited redshift range in which we can observe 
(z < 0.112). Nevertheless, iVciv does not appear to be 
well-correlated with either WHIM tracers (O VI, N V) 
nor lower ions (C III, H I) and thus may arise from both 
photoionization and shocks. We conclude that C IV is 
not a conclusive WHIM tracer, though it is a good probe 
of metal enrichment. Longer wavelength observations of 
AGN sight lines would improve the statistics of this ion. 

We made numerous detections of both C III and Si III 
in IGM absorbers (39 and 53, respectively). Column 
densities of both ions show reasonable correlation with 
A'hi, and the dM /dz distributions are similar, leading us 
to classify both ions as predominantly photoionized. 

Detections of both Si IV and Fe III are not as plenti- 
ful as the other low ions, and conclusions are harder to 
reach. Still, the majority of the Si IV detections (12/20) 
also show Si III absorption, and A^siin and A^siiv are 
very well correlated. We compare our observed line ra- 
tios with a set of CLOUDY models and find them consis- 
tent with a photoionized IGM with \ogU = —1.6 ± 0.7. 
Our other adjacent ion pairs (CIII/IV) appear together 
in only eight absorbers and their column densities are 
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not well-correlated. However, they are consistent with 
models with logU = —1.4 ± 0.5. Based on these diag- 
nostics, we conclude that typical IGM clouds have an 
ionization parameter of log ?7 ~ —1.5 ±0.5. Photoioniza- 
tion models with U ten times higher would be required 
to explain the observed A^nv /Nqvi and A'^civ / ^ovi ra- 
tios, and this bolsters our interpretation of those species 
as collisionally ionized. 

-^Foiii appears stronger than predicted by CLOUDY 
models by a factor of 100 or more when compared to 
either A^sim or iVciii- This may be indicative of a softer 
ionizing spectrum above ~ 3 ryd than was assumed in 
the models, but that explanation causes a set of other 
problems. More sophisticated modeling is required to 
explain these observations. 

4.1. Status of the Baryon Census 

With an unprecedented number of H I absorbers, we 
conclude that the Lya forest makes up 29 ± 4% of the 
local baryonic ma terial by mass. This is compati ble with 
previous surveys (iPenton. St ocke. fc Shulll[200l as weU 
as simulations ()Dave et all 11999. 20011). C III and Si III 
detection statistics suggest that 3-25% of the Lya forest 
material is enriched with metals at a detectable level at 
low z. 

At our most sensitive limit, logA^ovi ^ 13.0, we can 
further account for 8.6 ± 0.7% of the baryons residing 
in the WHIM at temperatures around T = 10^-^ K via 
O VI observations. Our N V detections suggest a smaller 
WHIM fraction, consistent with nitrogen's lower abun- 
dance and lower range of temperatures over which it is 
is abundant. However, the d/V /dz distributions of both 
WHIM tracers are fitted with a power law of index /3 « 2, 
implying that weak absorbers contribute as much to the 
mass fraction as strong absorbers. More sensitive FUV 
surveys will likely increase the WHIM fraction. 

Additionally, there is the issue of "Lya bias" . In Pa- 
per 1, we examined only strong Lya systems (WLya > 
80 mA), yet A^ovi was shown to have little or no cor- 
relation with A^Hi- In this work, we set out to minimize 
this bias by including all detected H I systems, no matter 
how weak. Inclusion of weaker Lya lines in the survey 
boosts O VI detection statistics and all related quantities 

{dAf /dz, ^^Iq^) by 20% compared with the previous 
survey. However, even this survey results in some Lya 
bias since H I detections are required at some level, and 
the finite S/N of the data limits this to Wr.vrv ^,15 mA 
m even the best data. The results of iTripp et SI ()2008D 
suggest an additional ^ 20% increase if O VI systems are 
measured with no corresponding H I detection. Combin- 
ing these two effects, we suggest that the Lya bias on 
WHIM surveys is a real effect at the level of 20-40%. 
Thus we suspect the true value of rJwHiM as detected 



via O VI absorption is probably closer to 10% of fif,. 

Even reliable WHIM tracers such as O VI and N V can 
only detect material that has been significantly enriched 
above primordial metallicity. Thus, there is a metallicity 
bias in this survey. If we take the very rough 3%-25% 
value from the Si III and other low ion estimates dis- 
cussed above as a typical fraction of the IGM enriched 
to detectable levels, that may boost the true fraction of 
material at T = 10^-^ K by 3-25 times (~ 30% - 250% 
of the baryons)! Clearly the latter value is unphysical, 
but fiwHiM ~ 30% is plausible and consistent with some 
simulations. 

But what about the balance of the baryon census? 
If condensed objects (stars, ga laxies, etc) account fo r 
~ 7% of cosmological baryons (jSalucci fc Persicl ll999^. 
then ~ 50% of the local baryons remain unaccounted for. 
Cooling of diffuse gas is fastest at lower temperatures, so 
it is li kely that there is a vast hot reservoir of gas in the 
IGM (|Cen fc Ostrikei<19 99: ShuU 2003). This hot WHIM 
(T > 10^ K) is inaccessible through FUV absorption 
lines, and the primary metal ions (O VII, O VIII) are only 
detected in the X-ray. While O VII and O VIII absorp- 
tion has b een seen associated with the Galaxy and cer- 
tain AGN (jNicastro et al.ll2005t iFang. Canizares. fc Yaol 
[200l), there are no confirmed measurements of these lines 
in the IGM. This remains a challenge for sensitive X-ray 
spectrographs of the future. 

The scheduled August 2008 installation of the Cosmic 
Origins Spectrograph on HST will provide a great oppor- 
tunity to study the low-z IGM. With greater than ten 
times the sensitivity of STIS and a resolution compara- 
ble to that of FUSE, high resolution, high signal-to-noise 
observations of faint quasars will become feasible in very 
reasonable observing times. High throughput will en- 
able us to observe many fainter AGN targets, many with 
very large pathlengths (z < 0.45). This will increase the 
available total IGM path length by perhaps an order of 
magnitude or more. With this increase will come a pro- 
portional boost in detection statistics of WHIM tracers 
such as O VI and N V as well as Lya, Si III, C IV, and a 
host of other astrophysical interesting lines. 
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